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ABSTRACT 

A new computational method that uses polynomial equations and dynamical systems to evaluate logical propo- 
sitions is introduced and applied to Godel's incompleteness theorems. The truth value of a logical formula subject 
to a set of axioms is computed from the solution to the corresponding system of polynomial equations. A reference 
by a formula to its own provability is shown to be a recurrence relation, which can be either interpreted as such to 
generate a discrete dynamical system, or interpreted in a static way to create an additional simultaneous equation. In 
this framework the truth values of logical formulas and other polynomial objectives have complex data structures: sets 
of elementary values, or dynamical systems that generate sets of infinite sequences of such solution-value sets. Be- 
sides the routine result that a formula has a definite elementary value, these data structures encode several exceptions: 
formulas that are ambiguous, unsatisfiable, unsteady, or contingent. These exceptions represent several semantically 
different types of undecidability; none causes any fundamental problem for mathematics. It is simple to calculate 
that Godel's formula, which asserts that it cannot be proven, is exceptional in specific ways: interpreted statically, the 
formula defines an inconsistent system of equations (thus it is called unsatisfiable); interpreted dynamically, it defines 
a dynamical system that has a periodic orbit and no fixed point (thus it is called unsteady). These exceptions are 
not catastrophic failures of logic; they are accurate mathematical descriptions of Godel's self-referential construction. 
Godel's analysis does not reveal any essential incompleteness in formal reasoning systems, nor any barrier to proving 
the consistency of such systems by ordinary mathematical means. 
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1 INTRODUCTION 



In the last century, Kurt Godel's incompleteness theorems |15 | sent Shockwaves through the world of mathematical 
logic. The conventional wisdom is that Godel's theorems and his interpretations thereof are correct; the prevalent 
discussion concerns what these results mean for logic, mathematics, computer science, and philosophy |l26l l20l [161 . 
But as I shall demonstrate here, Godel's theorems are profoundly misleading and his interpretations were incorrect: 
his analysis was corrupted by the simplistic and flawed notions of truth value and proof that have troubled logic since 
antiquity, compounded by his misapplication of a static definition of consistency to a dynamical system. Exposing 
these errors reveals that reports of logic's demise have been greatly exaggerated; we may yet realize the rationalist 
ideals of Leibniz and complete the logicist and formalist programs of Frege, Russell, and Hilbert. 

There are two key principles here: first, that proof in formal reasoning systems is an exercise in solving systems 
of polynomial equations, yielding solutions that are sets of elementary truth values; and second, that certain self- 
referential formula definitions are recurrence relations that define discrete dynamical systems (and in turn infinite 
sequences of basic solution sets). It is a corollary to these principles that all of the various syntactic results from 
such calculations make semantic sense as the truth values of formulas, including solution sets that are empty or have 
multiple members, and dynamic solutions that change with each iteration and depend on initial conditions. 

These principles are familiar and uncontroversial in the contexts of elementary algebra and dynamical systems; 
they apply just as well when the basic mathematical objects are logical truth values instead of ordinary numbers. If 
you understand how to do arithmetic in different number systems, what it means to solve equations, and how to deal 
with recursive constructions like the Fibonacci sequence, then you can understand Godel's mistakes. You will find the 
powerful new paradigm of dynamic polynomial logic, which is a continuation of the pioneering 19th-century work of 
George Boole |4, 5 |. Dynamic polynomial logic is grounded in the intrinsic unity of logic and mathematics. 

Relative to classical logic, dynamic polynomial logic is paraconsistent, paracomplete, and modal. In this algebraic 
framework the misguided principle of explosion is corrected: inconsistent axioms are shown to prove nothing instead 
of everything. Moreover the principle of the excluded middle is clarified; in classical logic this idea is applied incor- 
rectly, reflecting confusion between arithmetical and algebraic systems. Dynamic polynomial logic computes precise 
solutions that can be interpreted as alethic and temporal modalities. 

1.1 Godel's Argument 

Godel considered formal reasoning systems as described in Whitehead and Russell's Principia Mathematica |[3T| , 
which was an epic attempt to formalize the whole of mathematics. Godel's basic argument was that every formal 
reasoning system powerful enough to describe logical formulas, proof, and natural numbers (like PM) must allow 
the construction of a special formula that is semantically correct but syntactically undecidable: true by metalevel 
consideration of its content, but impossible to prove or disprove by mathematical calculation within the formal system 
itself. This special formula, denoted both [R{q);q\ and 17Genr in Godel's paper, asserts that the formula itself cannot 
be proven within the system. Thus follows the apparent paradox, which Godel described in this way (as translated in 
lEol): 

From the remark that [R{q)\q] says about itself that it is not provable it follows at once that [R{q);q] is 
true, for [R{q);q] is indeed unprovable (being undecidable). Thus, the proposition that is undecidable in 
the system PM still was decided by metamathematical considerations. The precise analysis of this curious 
situation leads to surprising results concerning consistency proofs for formal systems, results that will be 
discussed in more detail in Section 4 (Theorem XI). 

Godel's Theorem VI states that there must exist (in a formal system like PM) a formula such as his special [R{q)\q] 
which can neither be proven nor disproven within the formal system that contains it. His claim that his special formula 
[R{q)\q] is semantically true is presented in the text of his paper but is not called out as a theorem. Godel's Theorem XI 
states that the existence of such an undecidable formula renders the consistency of the enclosing formal system itself 
an undecidable proposition. 
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1.2 The Truth (Value) Is Complicated 



The surprising result from my analysis is that Godel's special formula is neither semantically correct nor semantically 
incorrect; instead it is exceptional in a particular way, relative to the expectation that a formula should have a definite 
elementary value. Such exceptions, which appeared to Godel as 'undecidability,' are features not bugs in formal 
reasoning systems: it is appropriate that some logical formulas cannot be proven to be simply true or false, because 
they are in fact neither Evaluating Godel's special formula is analogous to taking the square root of a negative number: 
such a square root cannot be proven to be any real number, because it is not any real number More elaborate data 
structures are needed to describe the values of some logical formulas, just as complex numbers are needed to describe 
the values of some arithmetic formulas. 

We shall consider four different exceptions that can be raised by logical formulas: the static exceptions of unsatis- 
fiability and ambiguity and the dynamic exceptions of unsteadiness and contingency. These exceptions represent four 
distinct types of 'undecidability,' and each type has a particular semantic meaning and gives characteristic syntactic 
results. Ignorance of these exceptions has plagued logic since ancient times. One can already find in the riddles of 
Aristotle's adversary Eubulides, who originated the liar paradox in the 4th century B.C., perfect demonstrations of 
unsatisfiable and ambiguous logical propositions f24]. Yet the same exceptions occur in algebra with ordinary num- 
bers, and in this context they are well-understood and not at all controversial. With the appropriate data structures 
and algorithms, these exceptional results become no more problematic for logic than V— 1 is for algebra. These four 
exceptions are consequences of the two key principles discussed next. 



1.3 Logical Equations Give Sets of Elementary Solutions 

The first key principle is that the axioms in a formal system (including the definitions of formulas) constitute a system 
of simultaneous equations; hence the proper way to describe the truth value of a formula is to give the solution-value set 
to its defining system of equations — the possible elementary values of the formula when all the equations are satisfied. 
As in general algebra, a system of logical equations can have zero, one, or more solutions: thus the solution-value set 
can be empty, have one member, or have many members. 

Using binary logic with the elementary values true and false (T and F), there are four possible solution-value sets 
(hence four different truth values): the set {T}; the set {F}; the set {T, F}; and the empty set {}. The first two sets are 
the definite or 'unexceptional' results: if the feasible set of values for some formula is {T}, then that formula is neces- 
sarily true (equivalently, it is a theorem); if the feasible set is {F}, then the formula is necessarily false (its negation is 
a theorem). The last two sets are exceptional: if the solution-value set is {T, F}, then the formula is ambiguous; if the 
solution-value set is {}, then the formula is unsatisfiable (and the underlying axioms are inconsistent). This style of 
categorization extends easily to sets of elementary values with more than two values, including infinite sets (like the 
natural numbers) and uncountable sets (like the real numbers). In general we consider truth values in the power set of 
the set of elementary values. 

As Boole explained in his Laws of Thought |5 |, logical formulas with binary truth values can be translated into 
polynomial expressions: the coefficients 1 and represent the elementary values true and false; symbolic variables 
represent basic propositions, with each variable x subject to the constraint x^ ~xio ensure that its only feasible values 
are and 1; logical conjunction translates as multiplication; logical negation translates as the difference from 1; and 
logical disjunction translates as a certain combination of addition and multiplication. Boole presented a complete 
algorithm to translate any logical formula into a polynomial with ordinary integer coefficients. 

For example, using Boole's original method the logical formula x^y translates as the polynomial xy — x+\ and 
the logical formula y x translates as the polynomial xy — y-\-\. Using these translations, constraining each formula 
to be true (i.e. to equal the polynomial 1), and separately constraining each variable to be either or 1 yields the 
following system of polynomial equations: 

xy — x-\-\ = 1 
xy-y+l = 1 

2 \^) 
X — X 

2 

y ^ y 

The only solutions to these equations for the variables {x,y) are (0,0) and (1,1). You can see thatx and y have the same 
value in each solution, so you might expect the biconditional formula x o y to be a theorem given the axioms x-^ y 
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and y — ^ jc. Using Boole's original method, the logical formula x o y is translated into the polynomial 2xy — jc — y + 1 ; 
substituting either solution for (x^y) the value of this polynomial is 1. Thus the solution-value set for the polynomial 
2xy — X — y + 1, subject to the constraints in Equation[Tl is the set {!}. 

Translating back into logical notation, the solution-value set for the formula x-^y given the axioms x^y and y^x 
is {T}; therefore x-^yv& called a theorem relative to these axioms. However, neither the formula x nor the formula 
y has a definite value given these equations; each formula has the solution-value set {0, 1} in polynomial notation, 
which is {F, T} in logical notation. Therefore each formula x and y is called ambiguous given the axioms x — > y and 
y — > X. Note that neither x nor ^x is a theorem given these axioms; likewise neither y nor -^y is a theorem. This 
'undecidability' for the formulas x and y is a correct interpretation of the relevant equations, not a sign of pathological 
incompleteness in formal reasoning. 

Infeasible equations have empty solution sets and render unsatisfiable all formulas that are subject to them. For 
example the simultaneous equations: 

z=\-z, z^^z (2) 

have no solution; the first constraint is violated for both values z = and z = 1 that satisfy the second constraint. 
Therefore the formula z is unsatisfiable subject to these equations, as are the formulas 0, 1, and 1 — z; all share the 
empty solution-value set {}. The equation z = 1 — z is a polynomial translation of the logical axiom z—^z which 
states that the formula z is defined to be true exactly if it is not true. This is one way to model the liar paradox, and in 
this version the problem is no more paradoxical than the unsatisfiable equation 1=0. 

1.4 Self-Reference Gives Dynamical Systems 

The second key principle is that, if it is permitted for a formula definition (or any other axiom or equation) to refer to the 
solution generated by evaluating its own system of equations, then the value of every formula may gain a more complex 
data structure. Such a solution self-reference is a recurrence relation that can be interpreted in two distinct ways: in 
the dynamic interpretation the recurrence is used as such to define infinite sequences of solutions governed by a 
discrete dynamical system; and in the static interpretation the recurrence is used to generate an additional simultaneous 
equation. Both interpretations can be reasonable and useful, but it is important not to confuse them during analysis. 

To illustrate, let us consider the quadratic equation 2x^ + 3x + c = in which we define the coefficient c to be the 
number of real solutions (forx) to the equation in which it appears. In other words we have the following specification 
for a system of equations, using real- valued variables x, c € M: 

2x2 + 3x + c = (3) 
c := the number of real solutions for X to Equation[3] (4) 

In this specification the value of c depends on itself. To model this dependence, let us introduce an evolution function 
denoted F{c) that gives the number of real solutions to Equation |3] when c takes the value supplied as the function's 
argument. In the dynamic interpretation we take the specification to define a system that changes over time according 
to this evolution function F{c), with the state Ct+\ at the next time f + 1 given by the value F{ct) of the evolution 
function applied to the current state c,. Thus we have: 

2x,^ + 3x, + Q = 0, c,+ i^F{c,); f e {0, 1,2, . . .} (5) 

Alternatively, in the static interpretation we take the problem specification to mean that the input value c and the output 
value F{c) must agree simultaneously, as in: 

2x2 + 3x + c = 0, c^F{c) (6) 

We must now determine what this evolution function F(c) is and which particular values of c need to be considered 
in solving the equations. For this problem the quadratic formula serves both needs; a more general approach will be 
presented later. As you may recall, the quadratic formula states that the number of distinct real roots of an equation 
flx^ + bx + c = depends on its determinant b^ — Aac: there are no real solutions if the determinant is negative, one if 
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it is zero, and two if it is positive. Therefore, the possible values of c are {0, 1,2}; and using the quadratic-equation 
coefficients a = 2 and b = 3 from Equation[3] the evolution function F(c) must satisfy: 



F{c) 



0, 9-8c<0 

1, 9-8c = 

2, 9-8c>0 



(7) 



Therefore F(0) = 2 since 9-80>0;F(l)=2 since 9 - 8 • 1 > 0; and F(2) = since 9 - 8 • 2 < 0. Using polynomial 
interpolation, a closed-form function F(c) can be constructed that performs exactly these mappings 1-> 2, 1 1-> 2, and 
2h^0: 



F[c) 



(8) 



For the static interpretation, substituting the evolution function from Equation [8] into the system in Equation |6] and 
then making explicit the domain of each variable produces the following system of equations: 



2jc^ + 3x + c = 0, c = -c^ + C + 2; ,y e M, ce{0,l,2}c: 



(9) 



After rearranging the second equation it is evident that this system has no solution: —2 implies c — V2 but c is 
required to be 0, 1, or 2. Thus in the static interpretation, the formula x is unsatisfiable and the whole system of 
equations is inconsistent. 

For the dynamic interpretation, substituting the evolution function from Equation |8] into the system in Equation|5] 
and leads to the following system of equations, which includes the recurrence Ct+\ ^F{c,): 



2x^ + 3x, + c,=Q, Q+i -c,2 + Q + 2; e R, qG {0,1,2} cM, f e {0, 1,2, . . .} 



(10) 



In this interpretation we consider c, to be the state of the dynamical system at time t. Each initial state cq G {0, 1,2} 
generates an infinite sequence (co,ci, . . .) of states at successive times f: 



^ (0,2,0,2,...) 
^ (1,2,0,2,...) 



co = 
CO = 1 

co = 2 ^ (2,0,2,0,...) 



(11) 



In the state q = the main equation 2xj + 3xt + c, = specializes to 2xf + 3x( + = which has the two solutions 

= which has two solutions 
which has no real solutions, 
for X, again depending on the 



X, = — I and Xt = 0. Then given Ct — I the main equation becomes 2xj + 3x, + 1 
Xt = — I and Xt = —\- Finally, given Ct —2 the main equation becomes 2x} + 3x, + 2 = 
Therefore Equation [TO] also generates a collection of sequences of solution-value sets 
initial condition cq: 

3 



co=0 
Co = 1 
CO = 2 



({-io},{},{-i,o},{},...) 
({-i,-^},{},{-i,o},{},...) 
({},{-i,o},{},{-|,o},...) 

These collections of infinite sequences are governed by a discrete dynamical system, 
compactly as a graph in which each node indicates a state of the parameter c and each 
set for the formula x that is generated from assuming the state corresponding to the orij 




(12) 

This system can be displayed 
edge shows the solution-value 
;inating node: 



(13) 



This graph, constructed from Equation[TOl offers some explanation as to why the static system in Equation|9]is unsat- 
isfiable: the matching dynamical system has no fixed points, only a periodic orbit that never reaches a steady state. 
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You can read off the sequences of states or solution-value sets generated by this dynamical system by following the 
edges in its graph. 

In general many different patterns could be discerned in the behaviors of dynamical systems; but for the purpose 
of classifying formulas we consider only the number of fixed points. If a dynamical system has exactly one fixed 
point then that system is called stable; if it has no fixed points then it is called unstable; and if it has more than one 
fixed point then it is called contingent. In these terms, the dynamical system defined by Equation [TO] and graphed 
in Equation [T3] is unstable. Unstable dynamical systems correspond to unsatisfiable static equations, and contingent 
dynamical systems correspond to ambiguous static equations. 

1.5 Anticipating Godel's Error 

We shall see that Godel's special formula [R{q);q] behaves just like the self-referential quadratic equation above. 
Using X to represent Godel's formula, it turns out that his definition specifies the evolution function F{x) : 1—x. When 
Godel's reference to provability is interpreted statically as the constraint jc = F{x), his self-denying formula [R{q);q] 
becomes an unsatisfiable system of equations: 

x^l-x, xe{0,l} (14) 
And when interpreted dynamically as Xt+i F{xt), Godel's formula becomes the recurrence: 

xt+i-i^l-x,, xe{0, 1} (15) 

This specifies a simple dynamical system with one periodic orbit and no fixed points. This dynamical system generates 
an alternating sequence of values for its state Xt at successive times f , for each initial condition xq: 

" xo = ^ (0,1,0,1,...) 1 
xo = l ^ (1,0,1,0,...) J ^ 

By generating these results, dynamic polynomial analysis will show that Godel's 'formula [R{q);q] that is true if and 
only if is not provable' is exceptional in precisely the same way as is 'the quadratic equation 2x^ + 3x + c = that has 
exactly c solutions': interpreted as static systems of equations, neither specification can be satisfied; and interpreted as 
dynamical systems, both specifications lead to sequences that oscillate infinitely and never converge to a fixed value. 

These types of 'undecidability' are not syntactic aberrations; they are semantically appropriate descriptions of 
the mathematical objects specified by their respective self-referential definitions. In fact, these complex results only 
seem exceptional because of the misguided expectation that the specified mathematical objects should have simple 
elementary values. Imagine the confusion that would result if we were to speak of 'the Fibonacci /ormM/a' (expecting 
it to have a definite numeric value) or 'the Fibonacci number' (expecting there to be just one) instead of 'the Fibonacci 
sequence' and 'a Fibonacci number' Although Godel's self-denying/ormM/a may seem puzzling, Godel's oscillating 
sequence and Godel's inconsistent equation are rather less so. 

It remains to be demonstrated that systems of logical axioms are indeed equivalent to systems of polynomial 
equations, and that self-reference of the type that Godel described is accurately translated using evolution functions, 
recurrence relations, and dynamical systems. But before discussing the details of translating logic to algebra, let us 
review several useful kinds of results that can be calculated by general methods of algebra from systems of polynomial 
equations and from discrete dynamical systems. 

2 GENERAL ANALYSIS OF POLYNOMIAL EQUATIONS 

Polynomials are an important class of formulas in elementary and abstract algebra. In this section we introduce 
terminology and notation to describe systems of polynomial equations and their solutions (for the task of formula 
evaluation), and we consider ways to construct and to count the members of a polynomial ring that meet certain 
desirable criteria (for the task of formula discovery). We develop simple algorithms to perform these tasks by hand for 
small problems, and make reference to more general and efficient methods from computational algebraic geometry. 
For the moment we are just discussing polynomials in the context of general algebra, without any mention of their 
provenance in logic and axiomatic formal reasoning. 



7 



2.1 Formula Evaluation 

We consider systems of polynomial equations whose variables and coefficients take values from an algebraic field 
K such as the real number system R, the rational number system Q, or a finite field F^/ of order d (with d a prime 
number). Note that in a finite field F,/ we must use integer arithmetic modulo d for calculation; for example in the 
binary finite field F2 the sum 1 + 1^0 since 2 = (mod 2). ifTTI provides a good general reference for polynomials 
and algebraic geometry; ifTOl and lfT2l describe fields and other structures in abstract algebra; modular arithmetic is 
discussed in ifTTl . 

Definition 1 (System of Polynomial Equations) Given an algebraic field K and a vector n := {x\,X2t ■ ■ ^x„) of vari- 
ables, the polynomial ring over the variables x with coefficients in K is denoted K[x]. We consider a set A of simulta- 
neous equations: 

A := {qi =0^ q2=0, q,„=0} 

in which each polynomial qj is a member of the ring /r[x] and also constrained to equal zero. The solution set to this 
system of equations, denoted '^{A) or 'f{qi = 0, ^2 = 0, . . . , qm = 0), is the set of values of x for which these equality 
constraints are satisfied, assuming that each variable x,- takes a value in the field K: 

r(A) := {xe/T" : <?i(x)=0, <?2(x)=0, ...,^,„(x)=0} 

By its construction the solution set 'V{A) QK" must be a subset of the ajfine space K". 

In the extreme case that the equations in A are inconsistent then 'f{A) is the empty set; conversely if the equations 
are tautological (or if there are no equations) then 1^(A) is the entire set K". In algebraic geometry, the solution set 
'V{A) is called the afftne variety defined by the polynomials {qi,q2, ■ ■ ■ ,qm} used in the equations. Note that in the 
special case n = that there are no variables, we imagine that the affine space has one member () that is the unique 
zero-length tuple. Thus with n = the solution set '/(O = 0) { ()} has one member since the equation = holds; 
however "/^(l = 0) {} is the empty set since the equation 1=0 does not hold. Additionally, we take the polynomial 
ring K[] with no variables to be the same as the original field K. 

We next consider the evaluation of a polynomial function subject to the constraints in a system of polynomial 
equations. Recall that the power set 2' of any set s is the set of all possible subsets of s (including by necessity the 
empty set and the original set s itself). 

Definition 2 (Solution- Value Set) Consider a system of polynomial equations as described in Definition [7] along 
with an objective formula p S K[x]from the same polynomial ring. The solution-value set J5^(/?) of the objective p 
subject to the equations A is defined as the set of its feasible values when the equations in the system are satisfied: 

^Aip) ■■= {pix)eK -.xeriA)} 

in other words the image under p of the solution set Y{A). By its construction the solution-value set S^a{p) C K must 
be a subset of the set of elementary values in the original field K. Equivalently the solution-value set must be a member 
of the power set ofK: 

yA{p) e 2^ 

Here using K to stand for the set of elementary values in the field as well as the field itself. 

If the solution-value set S^a{p) for a polynomial p is a singleton {k\ that contains just one member k Cz K, then 
we describe p as necessarily k subject to A; this is the unexceptional result. If the solution-value set is empty, then we 
describe p as unsatisfiable subject to A; this means that the equations A in the system are inconsistent. Otherwise the 
solution-value set must have more than one member (including the special cases that it is infinite or uncountable) and 
we describe p as ambiguous subject to A. 

At the extremes the solution- value set ^(75) may be the empty set or the entire set of elementary values from the 
underlying algebraic structure K. 

Moving on, when using coefficients from the binary finite field F2 it is possible to simplify systems of polynomial 
equations in the following way. 
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Lemma 3 (Conjunction of Binary Constraints) Any system of polynomial equations with coefficients in the binary 
finite field ¥2 or with each polynomial constrained to the values and 1 can be simplified to a single equation. For 
any a^b € {0, 1} the product ab = I if and only if both factors a — I and b ~ \. From this it follows that any set A of 
equations: 

{ q\ (x) = 0, ^2(x) = 0, . . . , ^„,(x) = } 

with each qj G F2[x] or each qj{x) £ {0, l}/or all feasible values ofx, must have the same solution set as the single 
equation: 

(^l(x)+l)(^2(x) + l)---(^„,(x) + l) = 1 

In this case the set A of equations given in Definition\l\can be replaced with the singleton {q* — 0} or the equivalent 
{q* + 1 = 1} whose solitary member uses the conjunction polynomial q* derived from the original equations: 

q* := (^i(x) + l)(^2(x) + l)---(^„,(x) + l)-l 

The conjunction polynomial produces the same solution set as the original equations: 

riq* = o) = riq* + i = i) = r(<?i = 0, <?2 = 0, ...,<?,„ = 0) 

and consequently identical solution-value sets for any objective formula. 

Note that in the special case m ~ that there are no constraints we take the product of zero factors to be the multi- 
plicative identity. Hence from A = {} we derive the conjunction polynomial <7* = 1 — 1 which gives the tautological 
constraint = 0. 

Example 4 To illustrate these definitions, consider the following system of two simultaneous equations in the vari- 
ables X and y (thus n — 2,x\ — x, X2—y,x. = and m = 2): 

x{y+\) = (17) 
y{x+l) = (18) 

Expanding these polynomials yields the constraint set: 

A := {xy + x^O, xy + y^O} (19) 

Let us proceed to evaluate this system of polynomial equations using two different algebraic structures: the real 
numbers R and the binary finite field F2. 

First we assume that the variables x and y and their coefficients take real-number values. Let us use T^(A)k to 
denote the solution set to Equations[l2]and[T8]using the algebraic structure R. According to Definition[T] 

nAh := {(-y)eM^: ^^t^Io} ^^0) 

You can see from inspection that only two pairs {x,y) of real numbers satisfy the equations, namely (—1,-1) and 
(0,0). Thus we have the solution set: 

r(A)R {(-1,-1), (0,0)} (21) 

This variety 1^(A)k is a subset of the affine space of all pairs of real numbers. Given this solution set, it is 
straightforward to identify the solution-value sets for the objective formulas x, y, and x — y+ I by substituting the 
values from 1^(A)jj into each objective: 

^Aixh ^ {-1,0} (22) 
^a(3')k ^ {-1,0} (23) 
^a(x-3'+1)m ^ {1} (24) 
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Next if we use the binary finite field F2 as the algebraic structure for the variables x and y and their coefficients, 
then the polynomial system defined by Equations [TtI and [TS] has the solution set: 

nAh. - (25) 
^ {(0,0), (1,1)} 

Recall that in the finite field F2 the sum 1 + 1 using modular arithmetic. Based on this solution set ^(A)^^, the 
solution-value sets for the objectives x, y, and x—y+ \ are: 

^a(x)f2 ^ {0,1} (26) 
yA{yh, ^ {0,1} (27) 
yA{x-y+l)w, {1} (28) 

In either polynomial ring or F2[x,y] the formula x — y+l has only one feasible value given the constraints in 

Equations [17] and [18] it is necessarily 1 . However, given these same constraints each formula x and y is ambiguous in 
either polynomial ring: the solution-value sets {—1,0} and {0, 1 } each have two members. 



2.2 Formula Discovery 

We now consider the inverse problem: subject to a set of polynomial equations, instead of computing the solution- 
value set for some polynomial objective formula, we must find the set of polynomials that yield a given solution-value 
set. 

Definition 5 (Inverse- Value Set) Consider a system of polynomial equations as specified in Definition\l\ along with a 
query set s QK of elementary values from the algebraic field K. The inverse-value set ,5^^ ' (i) is the set of polynomials 
in the ring K\s\ whose solution-value set is exactly s: 

When working with coefficients in the finite field F2, replacing the original set A of equations with the set {q* = 0} 
using the conjunction polynomial from Lemma\3\does not change the inverse-value sets. In that case, for every possible 
query s C K: 

=^{<,*=0}(*) = =^{<,i=0,^2=0....,</,„=0}W 

If the set A of equations is inconsistent then according to Definition^the solution-value set {p)for every polynomial 
p G K\%\ must be the empty set. Hence when A is inconsistent the inverse-value set =-5^' ({})/'"" empty-set query 
yields the entire polynomial ring A'[x], and the inverse-value set y^^{s) with any non-empty query s yields the empty 
set of polynomials. In other words, given infeasible constraints, every polynomial has no feasible value and no 
polynomial has any feasible value. 

Considering a fixed set of equations A, the inverse-value-set function ' is the inverse image or preimage of the 
solution-value-set function ■ Thus if S^^ ' (s) {/" 1 , /?2 , ■ ■ ■ } then for every polynomial pi in this inverse-value set 
it holds that S'^AiPi) — s', and furthermore the members {/?,} of this inverse-value set are the only polynomials in the 
ring /r[x] that have the solution-value set s. In algebraic geometry, the inverse-value set ^^({0}) with the query {0} 
is called the ideal of the affine variety '^{A): the set of every polynomial that vanishes on 1^(A). The characteristics 
of ideals ifTTl lead to the following special case. 

Lemma 6 (The Ideal from a Single Equation) In the special case that the set A of equations contains a single equa- 
tion q ~ that is satisfiable, then every member of the inverse-value set =^{^'-o}({0}) must be the product of some 
polynomial p from K[x\ with the polynomial q from the equation; conversely every such product p x q must be a 
member of the inverse-value set: 

^{,io}({0}) = {pxq:peK[x]} 
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The products p x q are not necessarily distinct for different values of p. Recall from Definition \5\ that infeasible 
equations yield empty inverse-value sets for non-empty queries: thus '^{^'-()}({0}) ^ {} if the equation q = is 
infeasible. 

Since there are dedicated methods in algebraic geometry to compute the ideal generated by a set of polynomials 
(in my terminology the inverse-value set for the query {0}), it is useful to define inverse-value sets for queries other 
than {0} in terms of the ideal. 

Corollary 7 (Incrementing Inverse Polynomials) For any query {k} that contains a solitary value kdK, the inverse- 
value set ,5^^^ {{k}) can be calculated by adding that value k to every polynomial in the inverse-value set for the query 
{0}; 

■^A\{k)) = {p + k:pey^\m) 

In the case that the inverse-value set =5^' ({0}) was generated from a single constraint as in Lemma^ then the same 
incremented inverse-value set is also given by: 

^{,io}(«) = {p-^q + k-.pem} 

This corollary allows closed-form description of inverse -value sets for singleton queries. 

Continuing on, it is useful to appreciate that the set of polynomials over a finite number of variables with coeffi- 
cients in a finite field is itself finite. 

Lemma 8 (Counting Polynomials in Finite Fields) The polynomial ring Fj[x] over n variables x := {xi^X2: ■ ■ ■ 
with coefficients in a finite field of order d (with d a prime number) contains exactly 

d'" 

distinct polynomials, each of which can be expressed as a unique sum Y^jCjtj of coefficients and monomials. In this 
sum each coefficient cj is a member of the field and each monomial (power product) tj is a product of the variables 
xi through x„ in which each variable is raised to a power between and — 1. There are d" possible monomials. 

Since (by Fermat's Little Theorem) the identity a'' = a holds for every element a e of a finite field of order d (with 
d a prime number), exponents higher than c/ — 1 are not required in monomials that use F^/. For example, using the 
finite field of order 2 and two variables x and y, the resulting polynomial ring F2[ji;,3'] has a set of 2^ = 4 possible 
monomials: 

{xV,xy, A^ A"} {^,x,y,l} (29) 

In any polynomial in this ring, each of these 4 monomials must be assigned one of the 2 coefficients in {0, 1}. Thus 
there are 2^ — 16 possible polynomials in F2[x,3']; they are listed in Table[T] Alternatively, using coefficients from the 
finite field of order 3 there would be 3 — 19,683 possible polynomials in the corresponding ring F3[x,3']. 

Example 9 We can enumerate all sixteen polynomials in the ring F2[x,3'] with binary finite-field coefficients, and 
partition them into the four possible inverse-value sets (for the queries i = {}, i = {0}, .? = {!}, and s = {0, 1}) given 
the constraints A = {xy + x — 0, xy + y = 0} from Equations [IT] and [18] Table[T] shows these sixteen polynomials and 
their solution-value sets; Definition[TO]below explains how the table was created. From the table you can see that there 
are four polynomials whose only feasible value is zero: 

■^/'({O}) ^ {Pl, P7, PlU Pl3} ^3Q^ 

{0, x + y, xy + y, xy + x} 



There are also four polynomials whose only feasible value is one 

{P2, p&, pn, p 

{ 1, x + y+ 1, xy + y+ 1, xy + x+ 1 } 



=^ { P2, P8, pn, PU} 
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There are no polynomials whose solution-value set is empty (because the constraints are consistent): 



'^a\{}) ^ {} (32) 
The remaining eight polynomials share the solution-value set {0, 1} of both binary finite-field values: 

^a"'({0,1}) {P3, P4, P5, P6, P9, PlO, Pl5, Pl6} (33) 

Thus, using the constraints in Equations [TTl and [TSl every polynomial p, e F2[x,3'] is assigned to one of the four 
possible inverse-value sets for polynomials with coefficients in the binary finite field. 

Some of these inverse-value sets can be derived in another way. Since we are using polynomials with coefficients 
in the binary finite field F2, Lemma|3]states that the original constraints xy+x — Q and xy +y = in Equations [TtI and 
[T8]are equivalent to the single constraint: 

{xy + x+l){xy+y+l) = 1 (34) 
Thus the conjunction polynomial q* is given by: 

q* := {xy+x+l){xy + y+l)-l x+y (35) 

using modular arithmetic. Definition |5] assures us that for any query s the inverse-value set computed from the con- 
junction polynomial is the same as the original: 

^{;vlv.o}(^) = ^{;vUo,.^+,.o}(^) (36) 

Now, Lemma |6] says that every polynomial in the inverse-value set -5^{^^,,=o}({^}) "lust be the product of some 
polynomial p in F2[x,}'] with the conjunction polynomial q* = x + y just computed: 

^{x+,=o}({0}) = {px{x + y):pe¥2[x,y]} (37) 

And indeed, multiplying any polynomial in the ring F2 [x,y] by the sum x + y will yield one of the four polynomials in 
'^A^({^}) listed in Equation[30] For example we have the following products (keeping in mind that the finite field F2 
uses integer arithmetic modulo 2 in which x^ = x, 1 + 1=0, etc.): 

{y+l){x + y) +> xy + x (38) 

ixy)ix + y) ^ (39) 

{xy + y){x+y) ^ xy + y (40) 

{xy + x+y){x + y) ^ x+y (41) 

Furthermore, Corollary|2]states that for every polynomial p in the inverse-value set ({O})- the sum p+l must 

be a member of the related inverse-value set -5^{^|j,=o}({^})- ^'^^ ^^^^ illustrated by comparing Equations [30l 

andllU 



2.3 How to Solve It 

The calculations required for formula evaluation and formula discovery using systems of polynomial equations with 
coefficients in finite fields can be performed by an extension of the truth-table methods of |2r| and |32|. It is practical 
to carry out these calculations by hand with pen and paper for small problems. Alternatively, there are sophisticated 
algebraic geometry methods implemented in several widely-available computer algebra systems that accomplish the 
necessary calculations in a more efficient, robust, and scalable manner (using coefficients from finite fields as easily 
as rational or real coefficients). These computational methods are derived from the Grobner-basis algorithms for 
solving polynomial equations that were invented by Buchberger in the 1970s ||6l. Commands to perform the requisite 
calculations in a computer algebra system are included in Section|4] Table-based inference is explained presently. 
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For a system of m polynomial equations with n variables and coefficients in a finite field with d elements, naive 
table-based inference requires the construction of a table with d'^' rows and d" columns. The resulting d'' +" entries 
must be exhaustively enumerated. The following definition introduces the value worksheet data structure and an 
inference algorithm based on it. Each row of the value worksheet contains the same information as a traditional 
logical truth table, but in a flattened form. 

Definition 10 (Table-Based Inference witli Finite Fields) For a system of polynomial equations with coefficients in 
a finite field of some prime order d, the solution set, solution-value sets, and inverse-value sets described in 
Definitions\l]\2\ and\5\can be computed using the following algorithm. 

(a) Construct the value worksheet W(A).- 

(i) Make a table W(A) with a row for every polynomial pi G Fj[x] and a column for every point x^. G (F^/)"- 
Following Lemma^ this table W(A) will have d^' rows and d" columns (excepting the row and column 
labels). 

(ii) For every polynomial pi G ¥^i[x] and every point x^t G (F^/)", use the value Pii^k) of that polynomial 
evaluated at that point as the entry at row i and column k of the table. 

(b) Compute the solution set Y {A) from the constraint polynomials A: 

(i) If d — 2 and m ^ I, replace the set A of equations with the singleton {q* — 0} using the conjunction 
polynomial described in Lemma\3\(this step may be omitted for pedagogic purposes). 

(ii) For each constraint (qj = 0) G A find the matching polynomial pi ~ qj at row i in the table W(A). Ex- 
amining this row i of the table, mark as infeasible every column kfor which the value Pi{s.k) ot row i and 
column k is not zero. If there are no constraint polynomials then all columns remain unmarked. 

(iii) The unmarked columns identify the feasible points that constitute the solution set y{A). If all columns are 
marked as infeasible then the solution set is empty. 

(c) Compute the solution-value set S^A{Pi) for every possible polynomial objective pi G F^[x].- 

(i) Add a column labeled 5^A{Pi) to the table 'W{A) for solution-value sets. 

(ii) For each polynomial p, at row i, the solution-value set is the set of table entries at the unmarked columns. 
Thus the solution-value set for each polynomial pi is the set of every value pi{s.k) taken by that polynomial 
at a feasible point x,t indicated as an unmarked column k. 

(iii) If there are no unmarked columns then every solution-value set is empty. 

(d) Compute the inverse-value set ,y^^ (s) for every possible solution-value set query s C F^/.- 

(i) Examine the entries in the solution-value set column S^A{pi)- For every unique entry s in that column, 
select the rows i\, ij, etc. such that S^A{Pi\ ) = ^> '^AiPii) — ^' '^^'^ ^o on. The corresponding polynomials 
Pi^, pi^, etc. constitute the inverse-value set .y^\s). 

(ii) The inverse-value set for every query s C Fj that does not appear as an entry in the solution-value set 
column S^AiPi) empty. 

Example 11 Table [T] shows the value worksheet W(A) constructed according to Definition [TOl from the system of 
constraints A = {xy +x = 0, xy +y = 0} in Equation[T9l the polynomial ring F2 [x,y] with binary finite-field coefficients 
is used. There are two ways to mark infeasible points. Using the original constraint polynomials we note that xy-\-x 
appears as /513 and xy + y appears as pn in W(A). Examining row 11 we see that /?ii(0, 1) 7^ 0, therefore we mark 
the column for (0, 1) infeasible. In row 13 the entry /7i3(l,0) 7^ hence we mark the column for (1,0) infeasible. 
Alternatively, using the conjunction polynomial q* = x + y gives equivalent results. The polynomial jc + y appears 
at row 7 in W(A) and in this row both entries pi{Q, 1) and pt{1,Q) have nonzero values; hence the corresponding 
columns are marked infeasible. 
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Table 1 The value worksheet W(A) for the equations A = {.xy + x — 0, xy + y — 0} using the polynomial ring F2[x,y]. The 
points (0,1) and (1,0) are marked infeasible; the solution set 'y{A)f^ contains the remaining points (0,0) and (1,1). 

' Pijx^y) p,(0,0) PijO,!) P,(1,0) PijlA) ^fa+x=0.;o+y=0}(P/) 
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Each solution-value set ,S^A{pi) at row / is the set {/?,(0,0),/?,(l, 1)} of values taken by the polynomial at the fea- 
sible points (the unmarked columns in W(A)). For example in row 8 we have p%=x-\-y+l and both /?g (0, 0) = 1 and 
/?8(1, 1) = 1; thus the solution-valueset ^(x+y^- 1) {1}. The inverse-value sets are given by selected rows: for ex- 
ample rows 1, 7, 11, and 13 share the solution-value set {0} so the inverse-value set ^5^' ({0}) is { pi , py, /?n , }. 
The inverse-value sets for the queries {1} and {0, 1} are constructed in a similar way. However, since no polyno- 
mial has an empty solution-value set (the entry {} does not appear in the last column of the value worksheet), the 
inverse-value set of the empty set is itself empty: =5^' ({}) ^ {}. 

3 DYNAMICAL SYSTEMS FROM REFERENCES TO SOLUTIONS 

Now let us consider the complexity added to a system of equations when it is permitted to refer to the solution of 
that system of equations within the system itself. This type of self-reference is an elaborate recurrence relation, 
which can either be interpreted as such to define a discrete dynamical system, or interpreted in a static way to provide 
an additional simultaneous equation. Although much of the terminology and notation for 'dynamical systems' is 
relatively new (following a resurgence of interest in the 1970s, especially in nonlinear and chaotic dynamical systems), 
the mathematical treatment of recursion is quite old. For example, what we know as the Fibonacci sequence has been 
studied in various guises since at least the Middle Ages; and methods for: finite differences and difference equations 
have been developed since the work of Newton and then Taylor around the turn of the 18th century |3|. The modern 
treatment of dynamical systems lfT4l[T8l dates from Poincare's work at the end of the 19th century. It happens that 
the dynamical systems that we will encounter in the study of logic are quite simple: discrete time, finite phase-space, 
first-order, autonomous, and usually linear. 

3.1 Extended Systems of Polynomial Equations 

In order to develop a computable representation of solution self-reference, we introduce several new features to the 
systems of polynomial equations described in Section |2] parameters, equation templates, iterative assignments, and 
solution references. The idea is that each constraint may be specified as a template instead of as a simple polynomial 
equation; the parameters of the templates are allowed to refer to solution sets and solution-value sets from the systems 
of equations in which they reside. 
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Definition 12 (Parametric System of Poiynomial Equations) Consider a system of polynomial equations as speci- 
fied in Definition\l\ with variables x ;= {xi,X2, ■ ■ ■ ,Xn) and coefficients in an algebraic field K; the polynomials in the 
system are members of the ring K\\\. We introduce a tuple :— {9i,92t ■ ■ ,9d) of parameters with the requirement 
that each parameter Oj must take a value in some specified set Uj. Thus the set U of possible values for is given by 
the Cartesian product: 

For each variable Xi and each parameter dj we add a data-type constraint T to identify the appropriate set of possible 
values: 

Tv, : Xi e K 
Tg^: 9 J e Uj 

Instead of using a simple set A of polynomial equations as described in Definition\l\ an extended polynomial system 
is specified using a set A{&) of parametric constraint templates and a set D of assignment templates, in addition to 
the above data-type constraints. In the set A(&) :— {ai,(X2, ■ ■ ■ ,CXm} each constraint template a,- is a function of the 
parameters 0." 

a,: <?,(©) = 

such that for any instantiation of the parameter values the template a,- reduces to an ordinary polynomial equation 
with qi{&) G K\fi\. In the set D :— {5e^ ,8e^,. . . , dg^} each assignment template 5j specifies the value to be assigned 
to the corresponding parameter dj at every iteration, using some parameter-updating function Xj whose arguments 
may include the solution set 'y{A{&)) to the extended system of equations being defined, as well as any parameter 
including 9j itself (here the double left arrow denotes assignment): 

Se^: Oj <= A,( 01, 02,..., e^", 1^(^(0))) 

In particular a parameter-updating function Xj may use the solution-value set ■5^a{&) (p) for some objective formula p 
or the cardinality \S^a(&) if) I of such a solution-value set. Although every parameter-updating function Xj that assigns 
a value to a parameter dj must return a value in the designated set Uj of possible values for that parameter, the 
updating functions are not required to be polynomial. If an explicit updating function for any parameter 9j is omitted, 
then the identity function Xj{9) : 9 is used as a default. 

According to Definition [12] the self-referential quadratic system in Equations [3] and H] can be specified as the 
following constraint and assignment templates: 

Zx'. X e R 

Tc: c e {0,l,2}cR 

«! ; 2x^ + 3x + c = ^ ^ 

Sc- c ^ |^{ai}W| 

Here the variable x takes real values and the possible values of the single parameter c are the set U := {0, 1,2} of 
integers between and 2. This parameter c is assigned the cardinality of the solution-value set for the objective x 
(subject to the equation aO by the updating function X\ (c) specified by the assignment template 5c. 

For a different example, a system to generate Fibonacci-like sequences can be specified with two parameters and 
a pair of assignment templates, without the use of any conventional variables or any constraint templates: 



8h 



a e Z 

^ ^ I (43) 

a ^ b ^ ' 

b <= a + b 



We next develop an algorithm to make explicit the discrete dynamical system implied by a parametric system of 
polynomial equations. 
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Table 2 The state-transition worksheet for the dynamical system defined by the quadratic equation 2x^ + 3x + c = with 
c solutions, as modeled by the parametric templates in Equation 1421 Each value c, specifies a state of the system. The 
set A(c/) contains the instantiation of the constraint template «! at that state c/. The solution-value set i5^(c,) W gives 
the solutions for the objective x subject to the instantiated constraint. The successor state F(c,) is calculated by 
processing the assignment template Sc which asks the size of this solution-value set =5^(c,)(-'^)- 
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Definition 13 (Extracting the Evolution Function) A parametric system of polynomial equations as described in 
Definition |72] encodes a functional relationship between the value of the parameter and itself Here we extract 
this evolution function, denoted F{&), in two useful special cases. First, in the case that each parameter-updating 
function Xj is a simple algebraic function of the parameters in — without reference to the solution set Y{A(&)) — then 
the evolution function F is given by the tuple of parameter-updating functions: 

F{9i,92....,9e) : 

This is the usual situation with recurrence relations in general algebra. 

Second, in the case that the parameter domain U is countable and finite then we can proceed with hypothetico- 
deductive analysis to derive F by computing a mapping 0; 1— > 0'- for every parameter value 0, S U. To express the 
hypothesis that the parameters have some particular value 0,-, each parameter 9j is assigned a constant value 
9] from its domain Uj. Using these values the constraint templates are instantiated into an ordinary set A{&i) of 
polynomial equations: 

A{&i) := { ^1(0,) =0, ^2(00-0, ^,„(0,)=O} 

The solution set y(A{&i))for these instantiated equations is computed by the table-based algorithm in Definition \10\ or 
by a general algebraic geometry method as appropriate. Following this solution subroutine the assignment templates 
in the set D are processed: an updated value 9j is computed for every parameter via its parameter-updating function, 
using as arguments the solution-value set 1^(A(0,)) just computed along with the input values of the parameters: 

0; <= A,(0i,02,...,0/', ^(A(0,))) 
The desired evolution function F must map the hypothesized parameter values to their updated counterparts: 

(0i,02,...,0f) ^ (0i',0^,...,0;) 

The set of mappings 0/ &'ifor every parameter value 0; € U completely defines the evolution function F. When the 
parameter domain U is a subset of the field K used as polynomial coefficients, then the evolution function F can be 
expressed in closed fonn using polynomial interpolation. 

To illustrate the first case in Definition [T3] the evolution function for the Fibonacci-like parametric system in Equa- 
tion|43]is simply: 

F{a,b) : {b,a + b) (44) 

using the parameter-updating functions Xi (a) : b and X2{b) : a + b specified by the assignment templates 5a and 5/, in 
Equation [43] 

The requisite calculations for the second case in Definition [13] can be organized in a state-transition worksheet 
constructed as follows. In each row of the state-transition worksheet we record an index /, some value 0/ from the 
set U, the instantiation A(0,) of the constraint templates at that value, the relevant feature of the solution-value set 
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'f{A{&i)) given those instantiated equations, and the value f (0,) of the successor that was computed from y{A{&i)) 
and 0, by processing the assignment templates according to Definition [13] Each row of the state-transition worksheet 
gives a specific value F(0,) of the evolution function for the argument 0,-. The complete evolution function F(0) can 
be represented as a simple transition matrix or as a state-transition table; or using the method described in Lemma[T9l 
below, F{&) can be specified as a polynomial with coefficients in a finite field. 

Table|2]shows the state-transition worksheet for the parametric system in Equation|42] which uses a single param- 
eter c with domain U = {0, 1,2}. In this case the relevant feature of the solution set is the solution-value set ^A(ci){^) 
for the objective x. This worksheet constructs the following evolution function: 

F:{0,1,2}^{0,1,2} 

The function-development method of Lemma[T9]below enables us to construct a polynomial function that matches any 
transition function F :U ^ U for a finite phase space U (after assigning integers to identify the states if they are not 
already numeric); it happens that F{c) : —c^ + c + 2 agrees with the mappings in Equation |45] (using the polynomial 
ring K[c] with real coefficients). 

Definition 14 (Tlie Derived Dynamical System) A parametric system of polynomial equations as described in Def- 
inition\T2\in turn defines a discrete-time dynamical system with state whose phase space is the set U of parameter 
values and whose evolution function is F(&) from Definitional] The set {0, 1,2, . . .} of nonnegative integers is used 
as the domain for the evolution (time) parameter t. The state &t of the dynamical system at any time t is the result of 
the t-fold composition of the evolution function F applied to the initial state ®q: 

0, fM(0o) 

For example, given state 0o at time t —Q the state of the system at time t — 3 is the composed value F^^'^{&q) — 
F{F{F{&q))). Following the usual conventions for discrete dynamical systems, a sequence of successive states 0;, 
F{&t), F{F{&t)), ond so on constitutes an orbit l[14\\18]l . Furthermore a state &t* is a fixed point in the dynamical 
system exactly if it is its own successor: &,* ~ F{&t* )■ We categorize a dynamical system by the number of fixed points 
it has. Let us say that a dynamical system is steady // it has one fixed point; unsteady ;/ it has no fixed points; and 
contingent if it has more than one fixed point. 

Note that this categorization concerns the state of the dynamical system indicated by the parameters (0i , 02i • • • 7 (^e) 
rather than the solution-value set for any objective formula per se; it could happen in an unsteady system that the 
sequences of solution-value sets for some formulas are nonetheless monotonous. Note also that if a system of equations 
has no parameters, then the dynamical system extracted according to Definition[T4]will have the phase space U = {()} 
containing one state which is the empty tuple (). In this case the state-transition function F{) : () is the identity function 
and the solitary state must be a fixed point; the system is trivially steady. 

By the categorization scheme in Definition [14] the dynamical system derived from Equation [42] for the quadratic 
equation 2x^ + 3x + c — is unsteady; as the graph in Equation [T3] shows there is a periodic cycle and there are no 
fixed points. In contrast the Fibonacci-like system derived from Equation [43] is steady; as we shall see there is one 
fixed point at {a,b) — (0,0) even though the orbits through the other points do not converge. 

3.2 Static and Dynamic Interpretations 

A dynamical system from Definition [14] can be interpreted in two different ways to evaluate an objective formula. In 
the dynamic interpretation, the orbits in the dynamical system are used to generate infinite sequences of solution-value 
sets for the objective (a sequential report). In the static interpretation, the evolution function is used to provide an 
additional static constraint limiting attention to the fixed points; what is reported is the union of solution-value sets for 
the objective from these fixed points (a stationary report). Such static and dynamic interpretations of a self-referential 
system of equations reflect subtly different views; it is fine to try either or both for any given problem. 
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Definition 15 (Collected Sequences of Solution-Value Sets) Consider a dynamical system derived as in Definitions 
\1 3\ and U4\ from a parametric polynomial system of equations as described in Definition\T2\ For any given objective 
formula p G K[x\ or p d &, the dynamical system specified by U and F{&) encodes an infinite sequence of solution- 
value sets for p using each state Gq € U as an initial condition. Each infinite sequence .y of solution-value sets 
results from solving the systems of equations A(&t) instantiated from the constraint templates at successive states 

©O,01,©2,....- 

^A(0)(/'|0o) (^A(0o)(/'):=n(0,)(^)'^A(02)(p), •••) 

With reference to Definition\T3\ when the objective p refers to a parameter dj then the initial (hypothesized) value Oj 
for the current state should be used to evaluate that objective, instead of the updated parameter value 0j. 

The collection S^* of solution-value set sequences is defined as a relation between the states in U and the infinite 
sequences that arise from them. For every state ©,■ G U the collection contains a mapping to the sequence with the 
respective initial condition: 



Such a collection is usually written in square brackets as illustrated in the examples. If every solution-value set in 
a sequence has exactly one member (which is always the case when the objective is a parameter), then the braces 
around those sets may be omitted to simplify notation: thus the sequence ({A;o}, {^2}, ■ ■ •) could instead be 
written {kQ,k\,k2, . . .). 

For the parametric system in Equation|42] Equation [TT] in the introduction already illustrated the collection 
of solution-value-set sequences for the objective formula c (with the set brackets omitted since each solution-value 
set is a singleton). Likewise Equation \V2\ shows the collection '^^^1}^^^ °^ solution-set sequences for the objective x. 
Note that you can read the sequences in either collection from the graph in EquationfTSl 

Considering the Fibonacci example, the parametric system in Equation|43]specifies a dynamical system with phase 
space U = I?, state © = and evolution function F(a,b) : {b,a + b). Using this dynamical system the collection 

of solution-value-set sequences for the objective parameter a is given by: 



(fl,/7)o = (0,0) 

(fl,/7)o=(0,l) 

(fl,^)o = (2,l) 



<^ (0,0,0,0,0,0,0,... 
^ (0,1,1,2,3,5,8,... 
^ (2,1,3,4,7,11,18, 



(46) 



The initial condition {a,b)o — (0, 1) yields the familiar Fibonacci sequence. Again since each solution-value set is a 
singleton the braces around sequence elements have been omitted. 

The dynamical system can instead be interpreted in a static way to examine its fixed points. 

Definition 16 (Simultaneous Solution for Fixed Points) Consider a parametric system of polynomial equations as 
described in Definition \12\ and its extracted evolution function F{&) derived by Definition 1731 in the special case 
that the domain Uj for each parameter dj /n := (0i , 02j ■ ■ ■ : ^i) O- subset of the algebraic structure K used for 
the conventional variables and polynomial coefficients in the problem. In this special case the solutions from the fixed 
points in the dynamical system can be computed in an alternative way by solving an extended system of static equations 
in which each parameter Oj is treated as an indeterminate along with the conventional variables inx: = {xi,X2, . . . ,Xn). 
The extended static system includes the data-type constraints X and a polynomial equation transcribed from each each 
constraint template a and from each assignment template 5: 







e K 








«; : 


qi 


= 




0; 


= A;(0,r(A)) 
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Each parametric constraint polynomial qi is now treated as a member of the extended polynomial ring: 



K[x\,X2, . . . ,x„, 01, 02, • • • , 0^] 

Note that in this alternative formulation, fixed points that give empty solution-value sets are not reported. 

By Definition [T6]the parametric quadratic system in Equationl42lyields the following system of static equations: 



: 

a\ : 
5c: 



2x^ 



X e M. 

c e {0,l,2}cl 

-3x + c = 

c ^ -c^ + c + 2 



(47) 



using the polynomial interpolation of the extracted evolution function F{c) shown in Equation |45] This system of 
equations is infeasible; thus the solution-value set =^{T^,T,..ai.5j(^) evaluates to the empty set, as does the solution- 
value set for any other objective formula subject to these unsatisfiable constraints. 

For the Fibonacci-like system in Equation |43] transcribing the assignment templates into static equations yields 
the following set A of equations (using the extended polynomial ring Z[a, fe] with integer coefficients): 



T„ : 


a 


e 


z 


Tfc : 


b 


e 


z 


5a: 


a 




b 


5b: 


b 




a 



(48) 



The unique solution {a,b) ^ (0, 0) identifies the fixed point noted earlier. Consequently for the objective a the solution- 
value set {0}. Interpreting its familiar recurrence as a static constraint, it is appropriate to say that 
'the Fibonacci number' is zero. 



Example 17 Consider a new quadratic equation 



3bx - 



n 

2 



in which the coefficient b is defined to be the 



number of real solutions to the equation in which it appears. This problem is specified as the following parametric 
system, using the conventional variable x and the parameter b: 



5b 



rrX - 



b e {0,1, 2}c: 



'3bx + \^ = 



(49) 



By Definition[T4]these templates yield a dynamical system with the phase space f/ = {0, 1 , 2} and an evolution function 
F{b) that must satisfy the following criteria: 



F{0)=0, F{1)=0, F{2)=2 



(50) 



The polynomial F{b) : b^ — b (constructed by Lemma [T9]l meets these criteria for b € {0, 1,2}. The derived dynamical 
system, labeled with the solution-value sets discussed next, is displayed in this graph: 



{} 



{-11,-1} 



(51) 



There are fixed points at b = and b = 2 and no nonconvergent orbits; in the terminology of Definition [14] this 
dynamical system contingent. By Definition[T5]the dynamical system yields the following collection of solution-value 
sets for the objective x: 



bo = ^ ({},{},...) 
bo = l ^ ({},{},...) 
bo^2 ({-11,-1}, {-11,-1},. 



(52) 
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Reviewing the fixed points and their associated solution- value sets shown in the graph in Equation IST] the dynamic 
interpretation shows that there are two cases in which the equation + 3bx + ^ = has exactly b real roots: when 
b~Q and there are no solutions, and when b ~2 and there are two solutions (namely — 11 and — 1 ). 

Finally, using Definition[T6]and the polynomial interpolation F{b) : b^ — b of the evolution function, the templates 
in Equation|49]can be transcribed into a static set of equations using polynomials in the extended ring M[.)c,/7]: 



5h 



X 

b 

-3bx + ^■ = 



e {0,1, 2}c: 



b 



(53) 



The static polynomial system in Equationl53]has two solutions (x, fe) = (— 11 , 2) and (x, /?) = (— 1 , 2); these identify the 
fixed-point solutions. Consequently in this static formulation the solution-value set for x is given by ^{T,..Tj,a[,5j}(-f) ^ 
{—11,-1}. Note that this static interpretation does not reveal that b = Q gives a valid instantiation of the template 
equation a\ , since in that instantiation there is no real solution for x. 



4 TRANSLATION FROM LOGIC INTO ALGEBRA 



Let us now visit the foundations of logic and explore the relationship between logical reasoning with truth values and 
mathematical reasoning with ordinary numbers. What we call logic — ^reasoning with binary truth values, formulas, 
axioms, theorems, and proof — is nothing more than algebra presented in peculiar notation. This basic equivalence was 
asserted by Boole in his innovative mathematical treatment of logic [4, 5 1. Indeed, it seems apparent that in both logic 
and algebra we find common concepts of number, operation and formula; variable, function, equation, and solution to 
equations; recursion and infinite sequence. However, through a cascade of historical misunderstanding and obfuscation 
(much of it associated with Godel's incompleteness argument), the intrinsic unity of logic and mathematics has not 
achieved universal acceptance. Especially, recurrence relations and dynamical systems have not been recognized as 
such in the contexts of logic and set theory. 

Considering these things I have extended Boole's ideas, and the two key principles presented in the introduction 
of this essay, into several A rf/cZei of Algebraic Translation for Logic: 

I. Logical formula Polynomial formula 

II. Axiom ^ Polynomial equation 

III. Solution-value set Truth value (when self-reference is forbidden) 

IV. Polynomial with solution-value set { 1 } Theorem 
V. Reference to provability ^ Recurrence relation 

VI. Features of discrete dynamical system ^ Truth value (when self-reference is allowed) 

Here the wavy arrow ^ should be read 'translates as,' with the idea that these general principles subsume many 
specific functional mappings. By analogy the general principle: 

Roman number Hindu- Arabic number 



includes the specific mappings IIII 4 and IV ~-+ 4. Godel's formula [R{q);q], which denies its own provability, is 
analyzed during the presentation of Articles|V]and|Vl] For concreteness, the examples in this section are accompanied 
by commands to perform their calculations in the computer algebra system Mathematica fyS\. Carnielli has studied 
Boole's polynomial formulation of mathematical logic and uses related ideas in his Polynomial Ring Calculus 17, JJ. 
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Table 3 Translation from logical to polynomial notation, to map logical formulas p and q into polynomials with real or 
finite-field coefficients. The atomic formulas are represented by variables x := (xi ,X2, . . . ,x„). When using real coefficients 
each Xi is subject to the constraint x? =x,- to limit its possible values to {0, 1}. 



Logical Polynomial in M[x] Polynomial in F2 [x] Description 



T 


1 






1 






Trae 


F 
















False 


-^P 


1 - 


P 




14 


P 




Logical negation (NOT) 


pAq 


pq 






pq 






Conjunction (AND) 


p®q 




q- 


2pq 


P^ 


-q 




Exclusive disjunction (XOR) 


pVq 


P^ 


q- 


pq 


P4 


-?4 


pq 


Nonexclusive disjunction (OR) 


p^q 


1 - 


p^ 


pq 


14 


■/?4 


pq 


Material implication 


p^q 


1 - 


p- 


q + 2pq 


14 


P^ 


q 


Biconditional (XNOR) 


p\q 


1 - 


pq 




14 


pq 




Nonconjunction (nand) 


piq 


1- 


P- 


q + pq 


14 


p^ 


q+pq 


Nondisjunction (NOR) 



4.1 Logical Formulas as Polynomials 

Let us consider two methods to represent logical formulas as polynomial expressions: the original method created by 
Boole, and a revised method using finite fields and modular arithmetic. In either case we start with a logical formula in 
the propositional calculus built from some atomic formulas {xi,X2, ■ ■ ■ ,x„) and the usual unary and binary operations 
(negation, conjunction, disjunction, material impUcation, and so on). Any such formula has a particular truth table, 
and it is truth tables that are directly translated into polynomial functions; in fact, each translated polynomial can be 
considered a closed-form representation of a truth table. Many different logical formulas might share a given truth 
table; hence many different logical formulas might have the same polynomial translation. Translation is not limited to 
2-valued logic; in general any finite number of elementary truth values can be accommodated. 

Although 1 discuss polynomial translation, for Boole polynomials were not used to translate from some other sym- 
bolic notation for logical formulas; they were his only mathematical representation. Notably, Frege's Begrijfsschrift 
was published many years after Boole's lifetime, as were the subsequent works from Peano and Hilbert in which the 
ideography of contemporary logical notation was developed. 

The detailed translation methods described below have been used to generate a table of polynomials representing 
common logical formulas, including the elementary truth values, the unary negation of a formula, and the common 
binary operations applied to a pair of formulas. These translations, which are often called the Stone isomorphisms 
after 1281 , appear in Table [3] A practical way to translate a simple logical formula is to apply the substitutions 
listed in the table recursively, until all traditional logical operations have been converted into polynomial form. Using 
TableOit is possible to translate any formula from the propositional calculus into a polynomial with either real-number 
coefficients, or coefficients in the binary finite field F2. 

For example the logical formula y A (z ® w), which is the logical conjunction of the atomic proposition y with the 
exclusive disjunction (XOR) of the atomic propositions z and w, can be translated to a polynomial in the ring M[w,y,z] 
by first substituting the polynomial form of the inner disjunction, then substituting the outer conjunction, and finally 
simplifying the whole expression using standard algebra: 

yA(z®w) yA(z4-w-2zw) (54) 

(z + w — 2zw) (55) 
^ yz-\-yw — 2yzw (56) 

Besides translating the entire formula y A (z ® w) , each variable in it must be constrained to limit its possible values to 
and 1 . This is accomplished by the equations = w,y^ — y, and 7^ = z (which would be tautological if we were 
using constraints in F2). 
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4.1.1 Boole's Original Representation Scheme 



Boole's method to represent logical formulas as polynomial formulas was presented it in preliminary form in Mathe- 
matical Analysis of Logic and in more complete form in Laws of Thought [4, 5|. This method produces polynomials 
with integer coefficients that can be manipulated using ordinary arithmetic. Boole used the number 1 to represent 
the logical value true, the number to represent the logical value false, and symbolic variables such as x, y, and z 
to represent atomic logical propositions. Boole noted that if a variable x is limited to the values and 1, then the 
equation x^ — x must hold; he identified this equation as a 'special law.' Based on this equation and its rearranged 
form.x(l — x) = 0, Boole devised a clever method to develop polynomial functions from arbitrary truth tables in two- 
valued logic. Each developed polynomial has the property that it agrees in value with its logical predecessor for any 
combination of truth values of the atomic propositions. 

Boole's function-development method is illustrated here for functions of two variables, which is the most useful 
case since traditional logical notation uses unary and binary operations (negation, conjunction, material implication, 
and so on). Suppose that we have two real-valued variables x and y, each restricted to values in the set {0, 1}. We 
require a polynomial function p{x,y) that yields some specified value zi (also either or 1) when x — 1 and y — 1; 
likewise we require p{l,0) = zi, p{0, 1) = Z3, and p{Q,Q) — Z4, with each value Zi S {0, 1}. These required values can 
be arranged as illustrated on page 76 of [Si, perhaps the earliest specimen of a logical truth table: 



X 


y 


p{x,y) 


1 


1 


Zl 


1 





Z2 





1 


Z3 








Z4 



Boole proved that the requisite function p{x,y) can always be calculated as the following polynomial, using the re- 
quired values Zl through Z4 as coefficients: 

p{x,y) := zixy + Z2x{l -y) + Z3{1 -x)y + Z4{1 -x){l ~y) (58) 

Following Boole's special law we add the constraints x^ — x and y-^—yto require that each variable must be either 
or 1. Note that for any x £ {0, 1} and y E {0, 1} only one of the four terms Equation |58] attains a nonzero value. This 
feature, which holds for functions of any number of binary variables, is the basis of Boole's polynomial translation 
method. 

For example, the translation of the logical exclusive disjunction x©y uses its truth table: 



X 


y 


xQ)y 


1 


1 





1 





1 





1 


1 












which produces the following polynomial according to Equationl58l 

pix,y) := xil-y) + il-x)y (60) 
x + y-2xy (61) 

along with the constraints x^ — x and y^ — y. It can be verified by simple calculations that p{x,y) satisfies the stated 
criteria: 1) = 0; /9(1,0) 1; /^(0, 1) = 1; and p{Q,0) ~ 0. Note that in Boole's original method there is no cause to 
invoke unusual rules of arithmetic such as 1 + 1 = 1. Boole's arithmetic was emphatically the standard fare, not what 
we now call 'Boolean algebra' . However, as we shall soon see it is helpful to use instead of real numbers coefficients 
in the finite field F2 and thus adopt integer arithmetic modulo 2 (in which 1 + 1 =0). 

Table[3]shows the polynomial representations of several logical formulas using Boole's original method, alongside 
their conventional forms. In general, for any formula of the propositional calculus in a logical system with atomic 
formulas x := (xi,X2, . . . ,x„), Boole's representation method yields a polynomial in the ring R[x] with real-number 
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coefficients. The coefficients are more specifically always integers; hence the Boolean polynomials can also be con- 
sidered members of the ring Z[x] with integer coefficients or the ring Q[x] with rational coefficients as convenient. 

Note that there are some polynomials in the ring Z[x] (hence also Q[x] and M.[x]) that do not correspond to any 
well-formed logical formulas at all using Boole's representation; Boole recognized this fact and called such polynomi- 
als 'not interpre table.' Specifically, in Boole's original scheme it is not the case that polynomial addition corresponds 
directly to either kind of logical disjunction (exclusive or nonexclusive): the polynomial x + y is neither the trans- 
lation of x(By (using XOR) nor the translation of xVy (using OR). The polynomial difference x — y is similarly not 
interpretable as a logical formula in Boole's original translation scheme. 

There were some minor flaws in Boole's presentation of his polynomial representation method. Although the 
general algorithm presented in Chapter V of his Laws of Thought is correct, Boole did not always apply his own 
algorithm correctly. For example on page 95 of |31| he translated y A{z®w) as yz + yw, having omitted the last term 
—2yzw. It is also confusing that Boole used the signs for addition and subtraction in two different senses (sometimes 
as the elementary arithmetic operators and sometimes to signify set union and set difference); either usage is fine but 
it takes careful reading to disambiguate the overloading. Finally, Boole did not have very robust methods to solve the 
multivariate polynomial equations that he had formulated. 

4.1.2 Revised Translation Using Linear Algebra and Finite Fields 

Instead of using an ordinary number system such as the real numbers M, logical formulas can be translated into 
polynomials with coefficients and variable values in a finite field F^/ of (prime) order d. At the same time the range 
of acceptable input formulas can be widened from those using binary logic to those using multivalued logic (with 
d ^2 elementary truth values). There are three main benefits to using coefficients in a finite field instead of the 
real numbers M (or the integers, rationals, or complex numbers). First, the polynomial ring ¥ci[x] over a finite set 
X := {xi,X2, ■ ■ ■ ,x„) of variables is itself countable and finite. Having the set of possible polynomials thus Umited 
allows a simple tabular approach to solving systems of equations. In contrast there are infinitely many polynomials 
in any ring M[x] (as well as Z[x] etc.). Second, polynomial translations of logical formulas are a bit simpler and more 
intuitive using finite-field coefficients; for example the operation of addition in a polynomial ring F2 [x] maps directly to 
exclusive disjunction in 2-valued logic. Third, every polynomial in a ring Ff/[x] corresponds to a well-formed formula 
in t/-valued logic; there are no longer any uninterpretable polynomials. 

First let us generalize the notion of a logical truth table into a finite-integer function. 

Definition 18 (Finite-Integer Function) A finite-integer function T : {TLdT ^ "^d of order d and arity n maps from 
the n-dimensional affine space (Z^/)" to the set Z^, where TLj^ := {0, 1,2, . . . ,ii — 1} is the ring of integers modulo d 
(with d ^ 2). Such a function T can be visualized as a table in which each row i describes the corresponding mapping 
{(^i, 1 ) '^1,2 J ■ ■ ■ J fli.n ) ^ Zh such that T (xi ,X2, ■ ■ ■ ,x„) — Zi when each xj = ai j. 



i 


XI 


X2 


Xfi 




1 


ai,i 


a\,2 ■ 


ai.„ 


Z\ 


2 


fl2,l 


a2,2 ■ 


a2.n 


Z2 


d" 


ad". I 


ad". 2 ■ 




Zd" 



Lemma 19 (Polynomial Encoding by Linear Algebra) Consider a finite-integer function T : (Z^)" "Lii of arity 
n and order d as described in Definition |7S] with the additional restriction that d is prime. This function T can be 
represented in closed form as a polynomial p over the variables x :— {x\ ,X2,... ,x„), with coefficients that are either 
rational numbers or integers in the modular ring Z^/, such that p{x) = T{x) for every point x G (Z^)" (using modular 
arithmetic to evaluate p if its coefficients are in Zi^i, otherwise ordinary arithmetic). The developed polynomial p 
requires d" coefficients, which can be computed as the solution of a system ofd" linear equations in d" variables. 

Referring to Definition U 8\ let us designate a point Hi :— (ai^\,ai2, ■ ■ • ,«;,«) G (F^^)" as an index vector; we consider 
the list (ai , a2 , . . . , Sid" ) of all d" possible index vectors to be arranged in some ( arbitrary) order which is used through- 
out this computation. By Lemma\8\there are d" possible monomials in the polynomial ring ¥^[xi,X2, ■ ■ ■ ,x„]; each of 
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them can be generated by using an index vector to supply exponents for the variables in the tuple x := {x\ ,X2,. . . ,x„). 
From each index vector aj we generate the corresponding monomial tj as the following product: 

tj := {xiPHx2p-^---{x„p-" 

The vector t:— {ti,t2, ■ ■ ■ ,td") contains all the monomials for the polynomial ring. An index vector a; can also be used 
as coordinates ft? r some point Xi € {^df ^ cis given by: 

X; := {aij,ai2T ■ ■ ,ai,„) 

Combining these uses, we construct a matrix M of size d" x d" in which each row M,- is the monomial vector t evaluated 
at the point x; = a; defined by the index vector a,-; 

M; t(a,) 

Equivalently the element at each row i and each column j of the matrix M is the value of the monomial tj defined by 
the exponents aj evaluated at the point x,- defined by the coordinates a,- ( taking O" = 1 for this calculation): 

Each matrix element is roughly (a,)^J, where exponentiation works element-wise and the resulting products are them- 
selves multiplied together 

Next we introduce the column vector z := {zi,Z2T--,Zd") in which each Zi '■= T{xi) is the value of the finite- 
integer function T evaluated at the corresponding point. Finally we use the column vector c := (ci ,C2, . . . ,c^n)for the 
coefficients of the polynomial p which are to be determined. The system of linear equations Mc = z is a restatement 
of the original function T, amended to include the computed monomial values. The desired coefficients c of the new 
polynomial p are given by the solution c = ' z to this linear system. The resulting polynomial p = c t= (M^ ' z) • t 
is the scalar product of the solved coefficient vector and the monomials. The matrix inversion and scalar product can 
be computed using ordinary arithmetic, in which case the inferred polynomial p will be a member of the ring Q[x] 
with rational coefficients (thus also a member q/M[x] with real coefficients); alternatively, the requisite calculations 
can be performed using integer arithmetic modulo d in which case the inferred polynomial p will be a member of the 
polynomial ring Z^[x] with coefficients in the ring of integers modulo d (which is the finite field when d is prime). 

Lemma[T9l assumes that the matrix M constructed by its directions is invertible. 

Corollary 20 (Encoding Indeterminate Polynomials) The polynomial construction in Lemma\]9\works essentially 
unchanged if the values in the finite-integer function T are left indeterminate. Using the elements of the vector z : = 
(zi , Z2 7 ■ ■ ■ , Zrf" ) of function values as symbolic variables, the result p — (M^ ^ z) • t inferred by Lemma [79] is a member 
of the extended polynomial ring F^[xi,X2, . . . ,Xn',Z\_^Z2, ■ ■ ■ ,ZdA that includes these d" new variables. This allows a 
fully-parametric closed-form representation of any finite-integer function T of order d and arity n as described in 
Definition \18\ at the cost of introducing d" new variables. 

Example 21 Using 2-valued logic (thus d — 2) and the variables x and y (thus x = {x,y) and n — 2) leads to the 
following d" — 4 index vectors according to Lemma[T9l 

ai = (l,l), a2-(l,0), a3 = (0,l), a4 = (0,0) (63) 

Using these index vectors as exponents generates the following monomials: 

f 1 = x^y^ = xy, t2 = a; =x, fs = x^y^ =y, f 4 = x^y^ = 1 (64) 

and thus the monomial vector t = {xy,x,y, 1). The following Mathematica commands compute these index vectors and 
monomials: 

d = 2; xs = {x, y}; n = Length[xs]; (* variables *) 

as = Tuples [Reverse [Range [0 , d - 1]], n]; (* index vectors *) 

t = Table [Apply [Times , xs~as[[i]]], {i, 1, d~n}] ; (* monomials *) 
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Continuing on, we construct the 4x4 matrix M according to Lemma [19] 



/ 1 1 1 1 \ 



M = 



10 1 
11 
\ 1 / 



(65) 



For example the second row M2 = (0, 1,0, 1) is the value of the monomials t — {xy,x,y, 1) at the point {x,y) — (1,0) 
defined by the second index vector 32- In Mathematica the matrix M is computed by: 

M = Table [ReplaceAll [t , Table [xs [ [j] ] -> as [ [i] ] [ [j] ] , {j , 1, n}] ] , 
{i, 1, d-n}] ; 

Now we add a vector z = (0, 1 , 1 , 0) of truth values for the logical formula x®y arranged in the same order as the 
index vectors in Equation|63] These correspond to the truth table shown in Equation|59l The desired coefficients in the 
inferred polynomial are given by the linear system Mc = z. Using ordinary arithmetic, this system yields the solution 
coefficients c = (—2, 1,1,0). Combining this solution c with the generated monomials t, the inferred polynomial p in 
the ring M[jc,3'] is given by c • t; thus p{x,y) : x + y~ 2xy. In Mathematica this same solution is accomplished by: 

z = {0, 1, 1, 0}; 

c = LinearSolve [M, z] ; 

p = c . t; 

For completeness we add the constraints x^ — x and y^ —y lo limit the values of these variables appropriately. This 
first linear-algebra solution agrees with the result from Boole's original translation method. 

Alternatively, solving Mc = z using integer arithmetic modulo 2 yields the solution vector c = (0,1,1,0) and 
consequently the inferred polynomial p{x,y) : x + y in the ring F2[x,y]. No additional constraints are necessary in this 
case. In Mathematica this same modular-arithmetic solution is computed by: 

c2 = LinearSolve [M, z, Modulus -> d] ; 
p2 = PolynomialMod[c2 . t, d] ; 

Leaving the desired values (zi, 22,23,24) indeterminate and computing (M^'z) • t produces the following paramet- 
ric polynomial pi in the ring R[x,y,zi ,22,23,24]: 



In factored form this is exactly Boole's polynomial shown in Equation|58] Using modular arithmetic instead yields a 
similar parametric polynomial p2 in the ring F2[jc,3',2i ,22,23,24]: 



In Mathematica these parametric solutions are produced by: 
Clear [z]; (* remove old numeric values *) 

zs = Table [Subscript [z , i] , {i, 1, d~n}] ; (* list of variables *) 
pi = Expand [ (Inverse [M] . zs) . t] ; 

p2 = PolynomialMod[(Inverse [M, Modulus -> d] . zs) . t, d] ; 

Using Lemma [19] and modular arithmetic, the familiar values and operations of 2-valued logic can be rewritten 
as polynomial formulas with coefficients in the binary finite field F2; the results are summarized in Table [3] The 
basic values true (T) and false (F) map to the respective elementary values 1 and 0. Other logical functions of atomic 
formulas x andy are translated as follows. Logical conjunction (AND) translates directly as multiplication: for x Ay we 
substitute the product xxy (abbreviated xy 01 x-yin the usual algebraic fashion). Exclusive logical disjunction (XOR) 
translates directly as addition: for x y we substitute the sum x + y. Logical negation translates as incrementation 
by 1 (or equivalently the difference from 1): for the negation ^x of any formula x we substitute 1 + x (which is the 



X3'2i — xyz2 — xyz3 + X3'24 — X2i + X23 — y2i + 3^22 + 21 



(66) 



P2 



xyzi + xyz2 + xyzi + xyz\ + X2i + X23 + ^21 + ^22 + 21 



(67) 
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same as x — 1 or x + 1 or 1 — x using integer arithmetic modulo 2; but not the same as —x with a unary minus, which 
is just X itself). Nonexclusive disjunction (the usual OR) is a polynomial sum: for x Vy we substitute x + y + xy. For 
the material implication x y we substitute 1 + x + xy. The biconditional x ^ y (also written x = y or designated 
XNOR) is translated into the polynomial 1 + x + y. Finally, the nonconjunction x f y (alternative denial, NAND, Sheffer 
stroke, x\y) is translated into the polynomial 1 +X3', and the nondisjunction x \.y (joint denial, NOR, Pierce arrow, 
Quine dagger) is translated into the polynomial 1 + x + y + xy. 

There is already a BooleanTable command in Mathematica to generate truth tables for binary-logical functions, 
and a Boole function to convert the logical values True and False to the numbers and 1. The Tuples [] statement 
above was written so that the index vectors are generated in the same order that the BooleanTable function uses (the 
variables in xs also need to be in alphabetical order, as the Sort command would produce). With these caveats it is 
simple to convert logical expressions into polynomials within Mathematica. For example the following commands: 

z = Boole [BooleanTable [Implies [x,y]]] ; (* truth table {1,0,1,1> *) 
pi = (Inverse [M] . z) . t; 

p2 = PolynomialMod [ (Inverse [M, Modulus -> d] . z) . t, d] ; 

produce the values 1 ~x + xy for pi and 1 + x + xy for p2 which you can see listed in Table [3] as the translations for 
X —i'y using real and finite-field coefficients. 

4.1.3 Translation as a Function 

The translation from a logical formula (in the propositional calculus) to a polynomial can be considered a function; 
every logical formula maps to a unique polynomial. But in general there are infinitely many distinct well-formed 
logical formulas that map to any given polynomial in a ring Fc/[x] with finite-field coefficients; in other words the 
translation function is surjective but not injective. By comparison, translation into Z[x] (hence also Q[x] or K[x]) is 
neither injective nor surjective — some polynomials with integer (hence also rational or real) coefficients correspond to 
many well-formed logical formulas and others correspond to none. 

Definition 22 (Logical Preimage of a Polynomial Function) Let us designate the set of well-formed logical formu- 
las that map to a given polynomial as the logical preimage of that polynomial; these formulas constitute an equivalence 
class in the sense that they always have the same value as one another ( their truth tables are identical). 

For example the logical formulas x y and ^x V y and (x — > y) A (^x V y) are each translated into the polynomial 
1 + X + xy in the ring F2 [x, y] ; therefore they are members of its logical preimage. 

4.2 Equations from Axioms 

An axiom in a logical system has the same meaning as an equation in an algebraic system: each is an assertion made 
for the purpose of deducing what follows. Equations and axioms are constraints rather than commandments: it is 
generally possible to write systems of axioms or other equations that cannot be satisfied. The essential nature of an 
equation does not change because its formulas originated in logical instead of polynomial notation. 

Logical axioms translate directly as polynomial equations, when those axioms are simple formulas from the propo- 
sitional calculus (not involving references to provability, quantifiers, or indeterminate predicates over infinite domains; 
any of which complicate matters as addressed later). Using the turnstile h to mark an axiom as in Frege's Begriff- 
sschrift ifTSl . an axiom h ^ is the assertion that the included formula q is true; in Frege's terminology the whole 
axiom is a 'judgment' and the included formula is its 'content.' Let us say that the logical system of interest has a 
set A :— {h qi,\- q2,... <?m} of axioms, in which each content qi is a formula built from some atomic formulas 
X := (xi,X2, . . . ,x„) and the usual logical operations. For concreteness we assume 2-valued logic, though this is not 
required for algebraic analysis. 

Using Table [3] or the general method in Lemma[T9leach axiom h qj can be translated into a polynomial equation 
with either real-number coefficients or binary finite-field coefficients. Each content qi is first translated from logical to 
polynomial notation using K[x] or F2 [x] as desired. Then each judgment is translated to a polynomial constraint in the 
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obvious way: the axiom h qi becomes the equation qi= 1 . The set A of translated axiom-equations has a conjunction 
polynomial according to Lemma[3] 

q* '■= (q'l)(?2)---(?m)-l (68) 

with a slightly different form since each right-hand-side is 1 instead of 0. The constraint q* =0 using this conjunction 
polynomial corresponds to the assertion that the logical conjunction qi Aq2 A ■ ■ ■ Aq,„ of the axiom contents is true. 
Therefore using binary finite-field coefficients a set A of axioms translates to a single equation featuring its conjunction 
polynomial: 

A ^ {^*-0}, <?*eF2[x] (69) 

Example 23 Lewis Carroll [8 1 provided this succinct version of his delightful logic puzzle about the barbers Allen, 
Brown, and Carr: 

There are two Propositions, A and B. 

It is given that 

(1) If C is true, then, if A is true, B is not true; 

(2) If A is true, B is true. 

The question is, can C be true ? 

Using the variables a, b, and c for the respective propositions A (that Allen is out of the shop), B (that Brown is out), 
and C (that Carr is out), we model the problem as the following axioms: 

h c^ia^^b) (70) 
h a^b (71) 

Translating the content of each axiom into a polynomial in the ring F2[fl,/7,c] with binary finite-field coefficients 
according to Tableland Lemma [19] yields the following polynomial equations: 

\- c-^{a^^b) abc+l = l (72) 

ha^b ab + a+l = l (73) 

By Lemma[3]the conjunction polynomial for this pair of equations is given by the following expression (calculated by 
integer arithmetic modulo 2 since the polynomials use coefficients in F2): 

q* := {abc+l){ab + a + l)-l abc + ab + a (74) 

Thus the axioms A in Equations |70] and |7T| are translated into to the following set of a solitary polynomial equation: 

A {abc + ab + a = 0} (75) 

using polynomial coefficients in F2. 

4.3 Solution-Value Sets as Truth Values 

Translating axioms from logical to polynomial notation is only the first step; next it is necessary to solve these equa- 
tions and interpret their solution. The truth value of a logical formula p subject to a set A of axioms is given by the 
solution-value set S^a{p) of that formula subject to those axiom-equations. In order to compute the solution-value set 
as specified in Definition|2l it is first necessary to compute the solution set "/(A) to the equations A that is specified in 
Definition[T] using either a general method from computational algebraic geometry or the manual table-based method 
given in Definition[lOl Calculations can be performed with real or finite-field coefficients. 
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According to Definition |2] the solution-value set ,S^a{p) C is a subset of the set of elementary values in the 
algebraic field K that contains the value of the objective formula p; therefore the possible values of ,5^a{p) are the 
members of the power set 2^. Thus using the binary finite field F2 with elementary values {0,1} there are four 
possible solution- value sets: 

{{},{0},{1},{0,1}} (76) 

It happens that these are the only four solution-value sets encountered even when logical formulas have been translated 
into the polynomial ring M[x] with real coefficients, since the range of each polynomial function is still limited to {0, 1 } 
by its construction in Lemma[T9] Anyway each set of elementary values in Equation|76]makes sense as the truth value 
of a logical formula. As described in Definition |2] these solution-value sets allow logical formulas to be categorized 
into those that are necessarily true (/') = {!})' necessarily false {p) — {Q}), ambiguous {.S^a (/') = {Oj 1 })' and 
unsatisfiable (^aI/?) = ID- 
Note that an objective formula is unsatisfiable if and only if the whole set of axiom-equations is inconsistent; in 
that case every possible formula is unsatisfiable subject to those same equations. For example given the unsatisfiable 
constraint 0=1 there is no feasible value for any objective formula: =5^jro=i}(0) ^ {}, '^{o=i}(l) ^ {}' ■^{0=1} W =^ 
{}, and so on. In contrast to the principle of explosion that an inconsistent set of axioms allows any formula to be 
proven true {ex contmdictione sequitur quodlibef), in the algebraic paradigm an inconsistent set of axioms allows 
no formulas to be proven true (nor false, nor ambiguous). Moreover, such undecidability by reason of inconsistency 
is distinguished from the undecidability by reason of ambiguity that results from trivial, tautological, or otherwise 
inconclusive axioms. For example ^{0=0} (■^) ^ {0, 1} differs from ^{0=1} W =^ {} since each is undecidable by a 
different mechanism. 

Example 24 Returning to Carroll's barbershop, it was established in Example |23] that axioms about Allen, Brown, 
and Carr translate as the equation abc + a/? + a = using polynomials in the binary finite field F2. The solution set 
y{abc + 0/7 + = 0) is easily calculated using the appropriate modular arithmetic in Mathematica: 

In[4]:= Solve[{a b c + a b + a == 0} , {a, b, c}, Modulus -> 2] 

Out [4]= {{a -> 0, b -> 0, c -> 0}, {a -> 0, b -> 0, c -> 1}, {a -> 0, 
b -> 1, c -> 0}, {a -> 0, b -> 1, c -> 1}, {a -> 1, b -> 1, c -> 0» 

The result contains five solutions for {a,b,c): 

y{abc + ab + a = 0) ^ {(0,0,0), (0,0, 1), (0, 1,0), (0, 1, 1), (1, 1,0)} (77) 

Using the constraints in Equations|72]and|73]directly (instead of their conjunction polynomial) gives the same solution 
set: 

In[5] := Solve[{a b c + 1 == 1 , a b + a + 1 == 1>, {a, b, c}, 
Modulus -> 2] 

0ut[5]= {{a -> 0, b -> 0, c -> 0}, {a -> 0, b -> 0, c -> 1}, {a -> 0, 
b -> 1, c -> 0}, {a -> 0, b -> 1, c -> 1}, {a -> 1, b -> 1, c -> 0}} 

You can verify that the logical back-translation of each solution satisfies the axioms in Carroll's problem stated in 
Equations |70l and [TT] For example with o = T, /> = T, and c = F as in the last solution {a,b,c) — (1,1,0): the 
material implication bis true; the inner material implication ^bis false; but since c is false the outer material 
implication c — > (o -^b) is nonetheless true. 

Carroll's problem statement requests the truth value of the logical formula c. It is evident from the solution set 
y{abc + 0/7 + = 0) shown in Equation|77]that both and 1 are feasible solutions for c. Thus we have: 

'^{.tc+ab+a=0}ic) ^ {0,1} (78) 

and it is neither the case that c is necessarily nor the case that c is necessarily 1 . The formula c is ambiguous — given 
Carroll's axioms the barber Carr could be either in or out of the shop. 
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Let us also calculate the truth value of the proposition abc. You can see by inspection of Equation [TT] that every 
solution in 'V{A) contains the value zero for at least one variable in {a,b,c). Thus the product abc must always be zero 
using these solutions; in other words the solution-value set for the objective formula abc is given by: 

y{abc+ah+a=Q}{abc) ^ {0} (79) 

This result indicates that the polynomial abc is necessarily 0, and therefore that the corresponding logical formula is 
necessarily false — given Carroll's axioms all three barbers cannot be out of the shop simultaneously. 

4.4 Discovering All Theorems 

In logic there is a special name for a formula that is necessarily true: it is a theorem. Therefore a logical formula 
/7 is a theorem given some set A of axioms exactly if its solution-value set ,S^a{p) has the value {1} (that is, the 
set containing exactly the number one). Furthermore the set of all theorems entailed by the set A of axioms (using 
some finite list x := (xi,X2, . . . ,ji:„) of propositional variables) is given by the inverse-value set ^5^' ({1}) described in 
Definition|5] Lemma|6l and Corollary|2l Every logical formula (whose propositional variables are x) that is a theorem 
given these axioms must translate as one of the polynomials in this inverse-value set (which is empty if the axioms are 
unsatisfiable). Using polynomials with finite-field coefficients, this algebraic formulation transforms the discovery of 
theorems from a search through an infinite set of logical formulas to a search among a finite set of polynomials (each 
of which corresponds to a unique truth table). 

It is left as a separate exercise to choose which member of the logical preimage of each polynomial theorem should 
be used to represent it; this is properly an optimization problem in the area of logic synthesis. For a polynomial with 
coefficients in the binary field F2, a reasonable default choice for its representative logical formula is the transliteration 
of the polynomial that maps multiplication back to logical conjunction (A, AND) and that maps addition back to logical 
exclusive disjunction (0, XOR). But note that it is not appropriate to use what is commonly called 'disjunctive normal 
form' since as Table[3]indicates, non-exclusive disjunction (V, OR) does not map directly to polynomial addition (using 
either real or finite-field coefficients). Note also that an exclusive disjunction xi 0X2® • • • ® ji:„ in 2-valued logic is true 
if and only if an odd number of its arguments are true; otherwise the disjunction is false. This corresponds to integer 
arithmetic modulo 2 which is appropriate for polynomials with coefficients in the binary finite field F2. 

The concept underlying theorems is not exclusive to logic; the more primitive notion is that a formula might have 
a definite solution value given some set of equations or other constraints. In general algebra the definite solution value 
zero is given special status. As mentioned after Definition|5] what is defined herein as the inverse- value set =5^7' ({0}) 
is in algebraic geometry called the ideal generated by the polynomials in the equations A. In this sense a theorem is 
almost ideal! Others have made the connection between logical theorems and polynomial ideals 81911221 . 

Example 25 Consider the axioms h x and h x — > y. Using Table [3]each axiom translates as an equation using polyno- 
mials with binary finite-field coefficients in the ring F2[x,3;]: 

h X X = 1 (80) 

hx-^y l+x+xy=l (81) 

By Lemma[3]these two equations yield the conjunction polynomial <7* = (x) ( 1 + x + xy) — 1 which evaluates to \+xy 
using modular arithmetic. Thus the axioms are translated as the set A = { 1 + xy 0} containing one polynomial 
equation. Using the table-based inference method in Definition [TO] or the Mathematica command: 

Solve [{1+ X y == 0}, {x, y}, Modulus -> 2] 

reveals that his set A of equations has the solution set y{A) = {(1,1)} for {x,y). Evaluating each of the 16 polynomials 
in the ring F2 [x, y] at this unique solution reveals 8 polynomials whose only feasible value is 1 (those whose listed value 
is 1 in the column labeled ^, (1, 1) in Table[T] namely p2, P3, ps, and so on). These 8 polynomials, which comprise the 
inverse-value set =^{^'+1=0} ({1})' the translations of all the theorems entailed by the stated axioms: 

1, y, X 1+x+y, xy, 1+y+xy, l+x + xy, x+y+xy (82) 



29 



Following Lemma|6]and Corollaryl?] this inverse-value set can also be described in closed form as: 

^{«+i=0}({l}) ^ {px{l+xy) + l:pe¥2[x,y]} 



(83) 



where modular arithmetic is to be used to evaluate the included polynomial expression. Each polynomial theorem has 
infinitely many formulas in its logical preimage; for example the last polynomial jc + y + xy corresponds to the logical 
formula {xAy)(Bx(By as well as to the logical formula V y. Using the matching formulas in Table [3] provides one 
choice of logical back-translation for each of the 8 theorems entailed by the axioms h x and h x — y: 



T, y, X, x-<r^y, xAy, y- 



■ X, 



■y, xVy 



(84) 



Included in this complete set of theorems is the sole formula y that a direct application of modus ponens would prove 
from the axioms h x and h x — > y. You can see that algebraic analysis provides more comprehensive results. 

Example 26 Returning once again to Carroll's barbershop, we can use algebraic analysis to discover the set of all 
theorems entailed by the axioms h c ^ (a ^ -^b) and h a ^ Z?. The conjunction polynomial q* = a + ab + abc in 
Equationl74lis a translation of these logical axioms. It follows from Definition|5] Lemma|6] and Corollary |7] that the 
set of theorem polynomials is given by: 



5^-1 

{a+ah+abc=Qi} 



({1}) 



{/? X (fl + ab + abc) + 1 : p G F2[fl,/7,c]} 



(85) 



By Lemma [8] there are 2^ = 256 distinct polynomials in the ring F2[a,^J,c]. Evaluating all of these polynomials 
with Mathematica reveals that there are exactly 8 theorems in the polynomial ring F2[a,fe,c] entailed by Carroll's 
barbershop axioms: 



{a+ab+ubc=Q} 



({1}) 



1 

1 + fl + fl/j + abc 

1 +flC 

1 + fl + fl/j + ac + abc 
1 -\-abc 
1+a + ab 
I +ac + abc 
I + a + ab + ac 



(86) 



The logical negation of each of these 8 theorems has the solution-value set {0}; thus each negated theorem is neces- 
sarily false. The remaining 240 polynomials in F2[fl,^7,c] are ambiguous, with the common solution-value set {0, 1}. 
Note that Carroll's original query c and its negation c + 1 are both members of this set -5^{7,+ai+„ic=o}('f^' ^ J') °^ 
ambiguous polynomials. 

For concreteness, the following Mathematica commands were used to generate the members of the inverse-value 

set '^[alab+ahc=o}i{^}) ^^own above: 



(* variables *) 
(* index vectors *) 
(* monomials *) 
(* all polys in ring *) 



d = 2; xs = {a, b, c}; n = Length [xs] ; 
as = Tuples [Reverse [Range [0 , d - 1]], n] ; 
t = Table [Apply [Times , xs"as[[i]]], {i, 1, d~n}] ; 
all = Tuples [Range [0, d - 1] , d"n] . t; 
(* construct inverse-value set from closed-form expression: *) 
theorems = Table [all [ [i] ] * (a + ab + abc) +1, {i, 1, Length [all] }] ; 
(* substitute to remove squares by Fermat ' s little theorem: *) 
theorems = Expand [theorems] /. {a~2 -> a, b~2 -> b, c"2 ->c }; 
theorems = DeleteDuplicates [PolynomialMod [theorems , 2]] (* modulo 2 *) 

The last command gives the output displayed above: 

{1, 1+a+ab+abc, 1+ac, 1+a+ab+ac+abc, 
1+abc, 1+a+ab, 1+ac+abc, 1+a+ab+ac} 
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4.5 Dynamical Systems from References to Provability 

A logical formula that refers to its own provability can be translated into algebraic form using a parametric system of 
polynomial equations as described in Definition [12] The parameters and parameter-updating functions introduced in 
the definition are expressly permitted to use features of solution sets and solution-value sets; this allows references to 
provabiUty to be modeled. From the specified parametric system of polynomial equations is extracted by Definition[T4l 
an evolution function F{&) that describes how the parameters change as a function of themselves. This evolution 
function can be used in one of two ways: either as a static constraint ~ F{@) that is added to the system of 
equations (which is then instantiated from parametric form into an ordinary system of equations; or as a recurrence 
relation 0,+i -4= F(0,) that is used to extend the static equations into a dynamical system. 

Example 27 Let us proceed now with the analysis of Godel's formula [R{q);q] that asserts its own unprovability. 
Using X to denote the formula and Bew(x) to denote the proposition that x is provable ('beweisbar' in German), we 
begin with the declaration x E {T, F} that the possible values of the formula x are the elementary logical values true 
and false. We define the formula x using the axiom: 

|-x = ^Bew(x) (87) 

Using Definition [12] and Table [3] the type declaration for x and the axiom in Equation [87] translate as the following 
parametric system of polynomial equations: 

T,: X e {0,1} CM 

5,: X l-(^{}(x) = {l}) 

For concreteness Equation [88] uses real-number coefficients; however for this problem the calculations would be 
identical using coefficients in F2 instead. 

Hypothetico-deductive analysis according to Definition [14] now reveals the evolution function F{x) for the para- 
metric system in Equation[88] Since the parameter x has the domain {0, 1}, we note that the function F must map from 
the set {0, 1} back to itself. We consider the two possible values of the parameter. In the case x — the parameter- 
updating function is instantiated as Ai(0) : 1 — (=5^{}(0) = {I})- Absent any constraints, the only feasible value of 
the objective formula is 0; in general for any constant k the solution-value set simply yields {k}. Hence 

the value of X\{Q) is 1 — ({0} = {1}), which is 1 since the false comparison statement yields the value 0. Similarly, 
in the case x = I, the parameter- updating function is instantiated as Ai(l) : 1 — (^{}(1) = {I})- The solution- value 
set o5^|}(l) ^ {1} gives Ai(l) 1 — ({1} = {1}) which evaluates to 0. These mappings 1— > 1 and 1 1— > define 
the evolution function F{x) encoded by the parametric system in Equation [88] It is evident from inspection that the 
interpolated polynomial function F{x) : I —x provides the appropriate values. 

There are static and dynamic ways to interpret Equation[88] In the static interpretation we replace the assignment 
template 5x with the static constraint x = F{x) using the extracted evolution function F{x) : 1 — x. In this static 
interpretation Godel's axiom in Equation[87]is translated via Equation[88]into this system of equations (the constraint 
could be expressed in polynomial form as —x using Boole's representation scheme): 

T.: X e {0,1}CR 
: X — 1 — X 

It is evident from inspection that Equation[89]has no solution; the solution-value set ^[T:^,s^-j (x) evaluates to {}. In fact 
simple substitution reveals that Equation[89]is equivalent to the constraint 1=0. Therefore, in the static interpretation, 
the truth value of Godel's formula [R{q);q] is that it is unsatisfiable. In other words, there is no such thing as a logical 
formula x that is true if and only if x is simultaneously not provable, just as there is no such thing as a quadratic 
equation 2x^ + 3x + c = that has exactly c real solutions for x. Each specification is internally inconsistent, and by 
algebraic analysis each inconsistency has been exposed as an infeasible system of equations. 

The dynamic interpretation of the parametric system in Equation[88]explains the mechanism of the inconsistency 
in Godel's formula. In this interpretation the parameter x, E {0, 1} gives the state of the dynamical system at each 
time t; the recurrence -4= 1 — .jc; from the evolution function F{x) : I —x creates a transition from the state x = 
to the state x = 1 and vice versa. The solution-value set for the objective x computed for each time t simply contains 
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the respective value of the state x,. Thus the dynamical system derived from Godel's formula is represented by the 
following graph: 

{0} 

®0® (90) 

{1} 

Following Definition [15] either initial value xq gives an alternating infinite sequence of solution-value sets for x at 
successive times t: 



xo = ^ ({0},{1},{0},{1},. 
xo=l ^ ({1},{0},{1},{0},. 



(91) 



This dynamical system derived from Equation[88lhas a periodic orbit and no fixed points; therefore it is unsteady in the 
terminology of Definition[l4] In this dynamic interpretation Godel's self-denying formula [R{q);q] specifies a discrete 
dynamical system that oscillates between the states of truth and falsity. The objective x describing the truth value of 
the formula oscillates between being necessarily false (having the solution-value set {0}) and being necessarily true 
(having the solution- value set {1}). 

Note that the set A of constraints is empty in the dynamic interpretation of Godel's formula; in the parametric 
formulation shown in Equation[88] there are no constraint templates a, only an assignment template 5^. In this respect 
Godel's formula is more like the Fibonacci recurrence than like the self-referential quadratic equations analyzed earlier. 

Example 28 It is illustrative to consider the complement of Godel's construction: a formula that asserts its own 
provability. Using the parameter y to represent the formula, we start with the declaration y e {T, F} and the axiom: 

hy= Be\N{y) (92) 

This type declaration and axiom are translated into parametric polynomial equations according to Definition [T2l 

Ty: y e {0,1} CR 

5y: y 4= {y^y{y) = {l}) ^^'^ 

Following Definition[T4l the identity function F{y) : y is extracted as the evolution function encoded by Equation|93] 
Incorporating the constraint y = F{y) the static interpretation of Equationl93]vields: 

j: y e {0,i}cR 
5y: y = y 

As the constraint 5y is tautological the solution-value set =5^{t,,,5,,}(3') evaluates to {0, 1}. Thus the formula that asserts 
its own provability is ambiguous in the static interpretation. 

In the dynamic interpretation the parametric system in Equation|93]specifies the following dynamical system based 
on the recurrence relation -4= y,: 

C® 0^ 

Either initial value yo gives a monotonous sequence of solution-value sets fory at successive times f: 



yo=0 ^ ({0},{0},{0},{0},. 
yo = l ^ ({1},{1},{1},{1},. 



(95) 



This dynamical system has two fixed points; hence it is contingent by Definition [14] In contrast to Godel's specifi- 
cation in Equation [87] that describes a formula that cannot be realized consistently, the complementary specification 
in Equation [92] describes a formula that can be realized in two different ways, either of which is consistent with the 
specification: in one case the formula is necessarily true (having the solution- value set { 1 }), and in the other case 
the formula is necessarily false (having the solution-value set {0}). The tautological vacuousness of self-affirmation 
mirrors the viciously-circular contradiction of self-denial. 
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Example 29 Let us consider a set of axioms that gives a slightly more complicated dynamical system than Godel's 
formula [R{q);q]. We now specify a logical formula z with the following properties: z is not provable; if zis ambiguous, 
then it is not true; and if z is not provable, the negation of z is not provable, and z is not ambiguous, then z is true. Thus 
in logical notation we have the declaration z € {T, F} and the following axioms: 



h ^Bew(z) 
h Viel(z) ^z 

h (^Bew(z) A^Bew(^z)A^Viel(z)) 



(96) 



Here the notation Viel(z) is introduced for the proposition that z is ambiguous ('vieldeutig' in German), in other words 
that its solution-value set ^a{z) — {0, 1}. 

We introduce the parameter 9 to represent the solution-value set for z- The domain of 9 is the power set 
{{}, {0}, {1}, {0, 1}} of the elementary values {0,1}. Following Definition [16] we translate the axioms in Equa- 
tion|96]into the following parametric polynomial system (here we have delayed polynomial translation of the material 
implication operators and represented negations in a slightly different way for convenience): 



as 



z 
9 

io = {i}) 

((0^{O,l})^(l-z)) 
((0^{l})(0^{O})(0^{O,l})^z) 

9 



e {0,l}cR 

e {{},{0},{1},{0,1}} 

= 

= 1 

= 1 

^{ai,a2,a3}(^) 



(97) 



This parametric system of polynomial equations specifies the following dynamical system (each node is labeled with 
a state 9, and each arc indicates the updated solution-value set calculated for z assuming the state corresponding to the 
originating node): 

^ {} {0} V_y 

Hypothetico-deductive analysis produces the following mappings for the evolution function F: 

{1}^{}, {0}k^{0,1}, {0,1} .^{0} (99) 
For example, given the hypothesis 9 — {} the constraint templates in Equation|97]are instantiated as: 
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({} = {!}) - 
(({}^{0,l})^(l-z)) = 1 
(({}^ {!})({} ^{0})({}^{0,l})^z) = 1 



which simplify to: 

(0)=0, ((0)^(l-z)) = l, ((l)(l)(l)^z) = l 
Translating the material implication operators according to Table [3] yields the polynomial constraints: 

= 0, 1 = 1, z=l 



(100) 



(101) 



(102) 



These constraints lead to the solution-value set ^{o=o,i=i,r=i}(z) =^ {1} for z, which is assigned as the new value of 
the parameter 9 by the updating function in template 5g. Thus is established the mapping {} iH- {1} for the evolution 
function F{9). The other hypotheses 9 = {I}, 9 — {0}, and 9 — {0, 1} are handled in a similar way. The results are 
the same using coefficients in M or in F2. 

For this problem there is no direct representation of the evolution function F{9) as a polynomial with coefficients 
in R or in F2, since the domain {{}. {0}, {1}, {0, 1}} for 9 is not a subset of either set of numbers. Nonetheless it is 
evident from the mappings in Equation|99]and the graph in Equation|98]that there are no fixed points in this dynamical 
system; there is no state of 9 that maps to itself. Hence the formula z defined by the axioms in Equation|96]is unsteady 
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in the dynamic interpretation and unsatisfiable in the static interpretation. Comparing the graphs in Equations [98] and 
|90]you can see that the dynamic behavior of this new formula z is a bit different from that of Godel's [R{q);q]. For 
example the new formula z produces two periodic orbits instead of one. Also, assuming that z is a theorem never leads 
to the conclusion that the negation of z is a theorem (as it does for Godel's formula); for the new formula z these states 
are in different, disconnected orbits. 

4.6 Quantifier Elimination 

It is not necessary to use quantifiers to examine Godel's incompleteness argument; however for the sake of com- 
pleteness let us address briefly how quantifiers should be handled in this algebraic formulation of logic. Existential 
and universal quantifiers add complexity to polynomial formulas and systems of constraints. In the field of compu- 
tational algebraic geometry, there is a large body of work on eliminating quantifiers from polynomial formulas with 
real-number coefficients ||29l l2l. Since Boole's translation method allows logical formulas to be converted into poly- 
nomials with ordinary numeric coefficients, such algebraic geometry methods for managing quantifiers are applicable 
to problems that were originally specified in traditional logical notation. Eliminating quantifiers from an equation may 
change the overall system of equations from a simple conjunction (all equations are to be satisfied simultaneously) to 
a disjunction (various subsets of equations must be solved separately). 

What is difficult about quantifiers is not the quantifiers themselves, but the indeterminate functions of infinite 
domains that often accompany them. Indeterminate functions of finite domains can be treated parametrically using 
Lemma[T9l However indeterminate predicates such as Man{x) and Mortal(.x;) assumed to have an infinite universe of 
discourse for x (as x e N) are problematic. Perhaps it would be possible to choose a 'quorum' (sufficiently-large but 
finite universe of discourse) based on the number and arity of predicates declared in a logical system. Anyway there is 
more work to be done to connect the treatment of quantifiers in computational algebraic geometry with the treatment 
of quantifiers in classical first-order logic. 

5 DISCUSSION 

I have argued two major propositions in this essay: that logic is an application of mathematics, more specifically 
solving systems of polynomial equations to yield sets of feasible elementary values; and that allowing equations to refer 
to their own solutions creates discrete dynamical systems. In this analysis I have employed Boole's groundbreaking 
algebraic formulation of logic, revised and extended with concepts from algebraic geometry and computer science. 
Within this framework of dynamic polynomial logic, Godel's special formula that asserts its own unprovability is 
easily expressed and is not at all paradoxical. Furthermore, Godel's formula is not so special after all: it is easy to 
modify ordinary quadratic equations to exhibit the same behavior, or to specify other logical formulas that encode 
more elaborate dynamical systems. Kurt Godel provided innovative but convoluted proofs of some normal properties 
of a simple recurrence relation, accompanied by spectacular misinterpretations. 

A familiar result from geometry may provide insight into the subtle cognitive error made by Godel and many 
classical logicians: what I shall call the Pythagorean fallacy. 

5.1 The Pythagorean Fallacy 

The Pythagorean theorem holds an important lesson about undecidability and incompleteness. We all know the famous 
equation + — relating the lengths of the sides and hypotenuse of a right triangle. Perhaps less well known is 
the fact that Pythagoras and his followers initially considered the only valid numbers to be ratios of natural numbers. 
To them, the solution for c with some values of a and b (such as a = 1 and b ^ I) was alogos (dXoyoc;): though we 
typically say 'irrational', a more literal translation would be 'not reasonable'. Of course, to the modern thinker it is 
not a number like \/2 that seems unreasonable; it is the expectation that an algebraic equation with natural-number 
coefficients must have a solution that is the ratio of two natural numbers. 

Let us first recognize the phenomenon of the anticipated-actual type mismatch: the expectation that the value 
of a mathematical expression should be one type of object when it is in fact another. It is an anticipated-actual 
type mismatch to expect the expression \J + b'^ to have a rational value for all natural-number arguments a and b. 
The anticipated-actual type mismatch is facilitated by expressions that use objects of one type to construct objects 
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of another type: for example whereas the sum or product of any integers must itself be an integer, the quotient of 
integers may not be. We can recognize that certain mechanisms of constructing expressions lead to certain types of 
mathematical values: for example if we are allowed the square root operation then we can make algebraic numbers 
from integers, and if we are allowed infinite series then we can make transcendental numbers from rational numbers. 
Notation and terminology may obscure the actual mechanisms in use. For example there is no square root sign in the 
equation + = c^, yet we must invoke this operation to solve for c. 

Let us call it the Pythagorean fallacy to encounter an anticipated-actual type mismatch but then to misplace the 
blame: to decide that the actual value, rather than the incorrect expectation, is at fault. This error is facilitated by for- 
mulas that have the expected type for some but not all arguments (in other words the expected type is a special case of 
the actual type): for example the value of v3^+4^ is indeed a rational number (5) but the value of + is not any 
rational number We may recognize certain circumscribed exceptions in actual values: for example whereas dividing 
one rational number by another generally yields a rational number, division by zero does not. It is not problematic to 
recognize such limited exceptions. But it is inappropriate to indulge the Pythagorean fallacy by declaring that some 
general mathematical method is fundamentally defective because of an anticipated-actual type mismatch, when the 
fault lies with an incorrect expectation. 

Ancient Greek mathematicians (and perhaps Pythagoras himself) eventually resolved their Pythagorean fallacy in 
an exemplary way: they broadened their concept of number to include the irrational values necessary for the general 
solution of algebraic equations like a^ + b-^ = c^. Our modern concept of number has grown to include many more 
'unreasonable' types of objects: not only algebraic numbers like V2, but also transcendental numbers like e and n 
and imaginary numbers like V"-T. We have come to understand that the mechanisms of geometric and algebraic 
construction ought to produce irrational numbers sometimes, even when the input values are natural numbers. We 
would find it absurd to accuse geometry or algebra of undecidability or incompleteness because there are expressions 
whose values are not ratios of natural numbers. 

5.2 Godel's Error 

Returning to the main results from ifTSl . what Godel's Theorem VI says is technically true but not detrimental to for- 
mal reasoning as he and many others have concluded. However it is emphatically not the case that Godel's special 
formula is semantically true; instead this formula that asserts its own unprovability is exceptional in a particular way. 
The types of 'undecidability' and 'inconsistency' that have been identified as exceptions — unsatisfiability, ambigu- 
ity, unsteadiness, and contingency — are semantically and syntactically appropriate features of systems of equations. 
Moreover, these exceptions can just as well be demonstrated in ordinary polynomial equations which had no origin in 
traditional logical notation. For example it is evident that the formula x has no decidable real-number value given the 
definition that x is a real number satisfying = — 1 ; and that the formula y has no definite real value when y is defined 
as a real number satisfying y^ = y; and that the formula z has no definite integer value when z is defined as an integer 
satisfying Zr+i 1 — Z;. Each of these formulas x, y, and z is exceptional in a particular way and should not have a 
definite elementary value. There is no need to divorce syntax from semantics. Undecidability is not a crime! 
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A EXECUTIVE SUMMARY 



In his famous 1931 paper Godel described a special logical formula that asserts of itself that it cannot be proven, in the 
context of a formal system PM that can express logical formulas, theorems, proof, and natural numbers. Godel claimed 
that this special formula must be semantically true but syntactically undecidable: impossible to prove or disprove by 
calculations within PM. Godel concluded that the existence of such a true-but-undecidable formula renders any formal 
system like PM essentially incomplete and incapable of proving even its own consistency. In this presentation I 
shall introduce a new method of analysis called computational algebraic logic that proves most of these conjectures 
wrong. Although Godel's special formula is indeed 'undecidable' in the sense that it has a value outside of the set 
{true, false}, such undecidabUity is a feature rather than a bug in formal reasoning: the formula is not semantically 
true. 

There are three main ideas to be discussed: exceptions, translation, and dynamical systems. First, mathematicians 
are unfailingly perplexed by formulas that give values outside of the sets they had originally expected; they give 
pejorative names to the unexpected exceptional values. Natural numbers are usually welcome. But others are called 
negative, fractional (broken!), irrational, or imaginary. If you knew nothing about complex numbers, you might say 
that \/— 1 is undecidable: the value of this formula cannot be proven to be any real number, even though its argument 
is real. Likewise ^/2, n, e, and so on seem exceptional and undecidable if you are expecting rational numbers. Godel's 
undecidable formula is just another unexpected object in this litany of un-natural numbers: the best solution is to find 
the right data structure to represent it, not to indict formal reasoning as incomplete and inconsistent. 

Second, translation from logical to algebraic notation helps to uncover the right data structure. You would not 
want to do your calculus homework or balance your checkbook in Roman numerals; likewise it is difficult to perform 
logical inference in traditional notation (our mash-up from Frege, Peano, Hilbert, and a few others). Boole already 
showed in 1854 that formulas in binary logic can be translated into ordinary polynomials, and that logical deduction 
can be performed by solving systems of polynomial equations. Boole's work is widely misrepresented; for example, 
he used ordinary arithmetic in which 1 + 1=2, not what is now called 'Boolean algebra'. Anyway, using Boole's 
original method, the logical axiom p translates as the equation p = I and the axiom p ^ q translates as the equation 
pq — p + \ = \, using real- valued variables p and q. The constraints p^ =p and q^ = q ensure that each variable must 
be either or 1. The truth value of the formula q, for which I introduce the notation S^{cj), is given by the set of 
solutions to these equations: 



This system of polynomial equations is easy enough to solve by inspection; the only solution is (;?,^) = (1,1) which 
gives the solution-value set S^(q) = {1}. So the only feasible value for q is 1; therefore, we call q a theorem. Voila, 
modus ponens by elementary algebra. What was missing from Boole's work was a general method to solve systems 
of multivariate polynomial equations. Luckily for us, Buchberger invented such a method in the 1960s; today his 
Grobner-basis algorithms are widely implemented in computer algebra systems. 

Next, what does Godel's special formula look like in algebraic notation? To define a formula x that is true if and 
only if X itself is not provable, we can use a real-valued variable x and the equation x = {.5^{x) ^ {1}) in which the 
idea 'x is not provable' is expressed by saying that the solution- value set S^{x) is not equal to the set {1} (we adopt 
the usual computer-programming convention that the test of inequaUty ^ returns the numeric value if its arguments 
are equal and 1 if they are not equal). This self-referential equation turns out to be a recurrence relation which is 
equivalent to the recurrence xt+i -^l—xt with each xt S {0, 1}. The mathematical object that is defined is a discrete 
dynamical system with two states (x = and x=\) and a periodic orbit that oscillates between them: 



The graph shows exactly what is exceptional about Godel's special formula: it defines system that has no fixed points 
(steady state). This dynamical system encodes for x an infinite sequence (. . . , {0}, {1}, {0}, {1}, . . .) of alternating 
solution- value sets, with a phase determined by what initial condition xq was assumed. 

Dynamical systems are the 'complex numbers' that capture the value of Godel's special formula and similar self- 
referential systems of equations. These dynamical systems are perfectly computable from self-referential formulas 
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by ordinary means and their existence does not jeopardize the consistency of the formal system within which they 
are derived (elementary algebra). I do not think that any mathematician or philosopher is troubled that the Fibonacci 
numbers are not all the same; we should understand that Godel's special formula is the same kind of construction 
and that it yields a similarly dynamic sequence of elementary values. It is no more correct to say that Godel's special 
formula is semantically true or false than it is to say that Fibonacci's formula is semantically zero, one, or two (or 
three, or five, or eight; you get the idea). 

Within this framework of algebraic logic, many other types of exceptions can be recognized through computation. 
To complement Godel's special formula we might introduce a new formula y that is true if and only if it is provable: 
y = {y{y) = {!})• This is equivalent to the recurrence relation yf+i <= y, with each y, e {0, 1}. What is defined is a 
discrete dynamical system with two fixed points: 



Such a dynamical system is contingent. It encodes for y two different infinite sequences of solution-value sets. From 
the initial condition yo = that y is not a theorem the sequence ({0} , {0} , {0} , . . .) confirms that y is never a theorem. 
And from the initial condition yo=^ that yisa theorem the sequence ({1},{1},{1},...) confirms that y is always a 
theorem. 

Rounding out the host of exceptional formulas are those whose solution-value sets are empty and those whose 
solution-value sets have multiple members. Using Boole's translation 1 —y for ^y, we might model the liar paradox 
with the equations z=l—z and = z for z € M; in this case ^(z) = {} and the formula z is unsatisfiable (just hke 
the equation = 1). We might add the complementary truth- teller problem using w = w and = w with w e M; in 
this case S^{w) = {0, 1} and the formula w is ambiguous. 

What Godel saw as undecidability does not show some foundational defect in formal reasoning any more than 
irrational or imaginary numbers invalidate arithmetic. The myth of incompleteness is dispelled when self-reference is 
unmasked as a recurrence relation that defines a dynamical system. Godel's incompleteness argument is commonly 
understood to prove that logic and mathematics are fundamentally incompatible, but in fact he demonstrated quite the 
opposite: for formal reasoning it is essential to recognize that problems in logic are problems in algebra, just as Boole 
had demonstrated. 
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B RUSSELL'S PARADOX: RECURSION REDUX 



The methods presented in this analysis of Godel's incompleteness theorems can be applied to another well-known 
problem in logic, Bertrand Russell's paradox (introduced in a 1902 letter to Frege that is reproduced in |30|). Russell 
asked whether a special set — the set of all sets that do not contain themselves — contains itself or not; either answer 
seems paradoxical. Although Russell's problem does not translate directly as a system of polynomial equations, it is 
nonetheless a system of mathematical constraints and the same insights apply. Russell's construction is unsatisfiable 
in the static interpretation and unsteady in the dynamic interpretation; it behaves exactly like Godel's special formula. 
Furthermore, a bit of reflection reveals that the specific issue with Russell's paradox is neither totality nor unrestricted 
comprehension but indefiniteness: his construction requires the set being defined to be used as a free variable in the 
predicate that defines it, which violates even the unrestricted axiom schema of comprehension. Regarding Russell's 
paradox, it is sets of humans not axioms of sets that have a problem with comprehension. 



B.l Algebraic Analysis 

Let us define S as a set of sets that includes at least some set r, and further define r as the set of all members of § that 
are not members of themselves. We desire the truth value of the formula r e r that asks whether the set r is a member of 
itself. The definitions of r and of § provide two constraints that constitute a generalized (beyond polynomial equations) 
system of constraints: 

r := {sGS:i^s} (103) 
r e § (104) 

In the special case that S is considered to be the set of all possible sets (assuming for the moment that such a construct 
is meaningful), this pair of constraints amounts to the classic formulation of Russell's paradox. But Equations [103] 
and II 04! are not limited to that special case; in particular § may denote a countable and finite set of sets and it is not 
important whether § contains itself. It is, however, essential to the problem that § contains r. 

In order to ascertain whether the set r thus defined contains itself, we use the solution-value set for the formula 
r e r subject to the constraints in Equations 1 1 031 and l 1 041 

(r e r) (105) 

Using the methods in Section [3] there are static and dynamic ways to interpret this solution-value set. In terms of 
Definition [12] we use the value of the set-membership expression r e r as the parameter Q for this problem. We note 
that the expression r G r takes a value in the Boolean finite field F? := {0, 1} according to whether the relation in it is 
false or true. Thus if the relation r € r is false we have (r e r) = (equivalently r ^ r) and 0=0; conversely if the 
relation is true we have (r G r) = 1 (equivalently r e r) and = 1. 

Using a hypothetico-deductive approach we consider these two possible cases = and = 1 for the parameter 
= (r G r) as illustrated in Table|4] In each instantiation A (0,) the conjectured value of the formula r G r is specified 
as an explicit constraint (since it cannot be substituted directly into an existing constraint as in the earlier examples). 
From either conjectured value we compute the other. That is, from the conjecture = 0i = indicating r ^ r the 
constructed set r includes itself and hence the computed value di = {r G r) ^ 1. Likewise from the conjecture 
= 02 = 1 indicating r E r the constructed set r does not include itself and so the computed value 02 = (r G r) 0. 
Thus from Equations 11031 and 1 1 041 we derive the state transitions i-^- 1 and 1 i-^- defining the following dynamical 
system: 




(106) 



with the states labeled r ^ r and r G r instead of = and = 1 for convenience. Since this dynamical system has no 
fixed points, we conclude in this dynamic interpretation that the formula r G r is unsteady subject to the constraints in 
Equation I 1031 and 1 1 041 Moreover, in the static interpretation the formula r G r is unsatisfiable given the constraints in 
Equation I 1031 and 1 1 041 since its conjectured and computed values never agree. 



38 



Table 4 The state-transition worksheet for Russell's set-inclusion query using the parameter := (r e r). 



i Oi 
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{/- = {ie§: 


: s ^ s}, r e §, r ^ r} 


{1} 


1 


2 1 


{r = {ie§: 


: i ^ i}, r e S, r e r} 


{0} 






Russell's special formula specifies exactly the same dynamical system as Godel's special formula: an oscillator 
with two states. This coincidence is not surprising: when the definition of a formula includes a recursive reference to 
a solution-set feature that has two possible states, the only way to make an unstable/unsatisfiable system is to have the 
conjecture of either state lead to the computation of the other Such instability is exactly what seems paradoxical about 
formulas like Godel's and Russell's (for some reason the contingent/ambiguous dual of each problem is not perceived 
as paradoxical). 

B.2 Russell's Free Variable Revealed 

In more conventional language, Russell's construction necessarily employs the set being defined as a free variable 
in the predicate that defines it — which renders that predicate indefinite. It is important to recognize such indefinite 
predicates, which Zermelo recognized as problematic in his axiomatization of set theory f34l|. In particular, when 
specifying the construction of a set by the action of a predicate upon the members of some universe of discourse, the 
members of that universe of discourse must be considered when evaluating the predicate for free variables. 

Lemma 30 (Revealing Free Variables in Set-Building) Consider the definition of a set y using some universe of 
discourse z and some predicate ^ {x): 

y := {x€z:0(x)} 

Every member of the set z that is a variable is in fact a free variable in the predicate 0(x). In particular, if the set y 
being defined is a member of the set z used in its own definition, then the instantiation <j) {y) of the predicate ( with y 
as the argument) must be evaluated during the construction ofy. It is evident that y is a free variable in this instance 
(piy). Leaving the universe of discourse anonymous as in the unrestricted definition y {x: 0(x)} does not change 
this property; if y is a possible value for x then the variable y must still be considered free in the predicate (j){x). 
Actually, it is not just the predicate <j){x) that is at issue; it is the entire right-hand-side expression {x € z '■ (j) (x)} in 
the definition. Even with a definition such as y := {x Cz z '■ T} whose trivial predicate T (the elementary value true) 
has no variables, the definiendum y must be considered a free variable in the definiens {x G z : T} if and only ify £ z; 
similarly, the definiendum y is a free variable in the definiens {x : T} of the unrestricted definition y : — {x : T}. 

It is important to recognize that Zermelo's |34| axiom schema of separation (Axiom III) contains two clauses: a 
separation clause and a definiteness clause. That axiom is repeated here: 

Whenever the propositional function €{x) is definite for all elements of a set M, M possesses a subset Mg 
containing as elements precisely those elements x of M for which (£(x) is true. 

The idea of a definite propositional function (addressed in LemmalSOll already entails the principle of separation, in the 
sense that unrestricted set comprehension always implies an indefinite predicate. For example, using the unrestricted 
definition y :— {x : (p{x)} the variable y must be free in the predicate (x) by the argument above; hence the predicate 
is indefinite. To be sure, an expression that includes free variables could still have a definite value if other constraints 
in the system limit the values of those variables; for the moment we assume the absence of such additional constraints. 

However, requiring separation without requiring definiteness does not close the loophole in naive set theory. Sim- 
ply introducing some independently-defined set-of-sets z to restrict set comprehension as in y := {x € z : 0(x)} is not 
sufficient to ensure definiteness; in the case y e z the variable y remains free in (j){x) and thus the predicate remains 
indefinite (at some point in the construction of y the predicate (j) (y) must be evaluated to determine whether to include y 
in itself). That is all to say that an unrestricted axiom schema of comprehension, retaining the definiteness clause from 
Zermelo's Axiom III but omitting its separation clause, would be sufficient to prevent Russell's paradox. Conversely a 
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restricted axiom schema of comprehension retaining the separation clause but omitting the definiteness clause would 
not prevent the paradox. The principle that a definition requires a definite expression in its definiendum is not at all 
peculiar to set theory! 

Returning to Equation |103l in the definition r := {5 e S : ^ s}, the variable r is free in the predicate s ^ s exactly 
if r is a member of §. Hence using this restricted definition it would violate the definiteness clause in Zermelo's 
Axiom III to consider r a set exactly if r e §. In the unrestricted definition r := {s : s ^ s}, the variable r is always 
free in the predicate s ^ s since r is considered to be a member of the universe of discourse from which the values of s 
are drawn. Hence both clauses of Zermelo's Axiom III would be violated by the unrestricted definition. Using either 
the restricted or the unrestricted formulation, at some point the predicate to be evaluated must become the expression 
r ^ r in which r is plainly a free variable; thus the predicate is indefinite. The precise issue with Russell's paradox is 
neither totality nor unrestricted comprehension; it is the attempt to define a set using an indefinite predicate. 

One can imagine a set theory in which Zermelo's definiteness clause is relaxed and thereby it is allowed to use 
the set being defined as a free variable in its own definition. In such a liberalized set theory it would be possible to 
use recurrence relations to define dynamical systems that give infinite sequences of sets, just as it is possible to use a 
recurrence relation to define the Fibonacci sequence in elementary algebra. Russell's construction would give one of 
these dynamical systems, characterized by a periodic orbit that oscillates between two states as shown in the graph in 
Equation ll06l The unrestricted 'set of all sets' (as fromy := {x : T}) would define a different dynamical system, as if 
by the recurrence relation yt+i U {yt} (this system would have no fixed points, but orbits that are neither periodic 
nor convergent). As in the earlier examples, these dynamical systems could instead be interpreted as static constraints 
(both Russell's special set and the 'set of all sets' y: = {x : T} would be unsatisfiable in this interpretation). 

The method of algebraic analysis detailed above illustrates exactly which exceptions can occur if a set y is allowed 
to be a free variable in its own definition (either from the explicit circumstance y E zoi from unrestricted comprehen- 
sion): unsteadiness and contingency. (Note that a self-referential set y could also have a steady, unique value and be 
unexceptional.) In this sense we gain the ability to discern more details about the various paradoxes in naive set theory, 
for example to clarify how 'the set of all sets that do not contain themselves' differs mathematically from 'the set of 
all sets' not otherwise specified (the constructions give structurally distinct dynamical systems, one with a periodic 
cycle and the other without). We can also discover similarities between exceptional constructions in set theory and 
exceptional constructions in the rest of logic and algebra, for example the congruence between Russell's 'set of all sets 
that do not contain themselves' and Godel's 'logical formula that asserts its own unprovability.' 

B.3 Functional Programming 

As a final thought on Russell's non-paradox, it is interesting to me as a computer programmer that the passage in 
Frege's Begriffsschrift that inspired Russell to create his set-inclusion problem also has a computer-programming 
interpretation. This is how Russell introduced his problem [23 |: 

There is just one point where I have encountered a difficulty. You state (p. 17) that a function, too, can act 
as the indeterminate element. This I formerly believed, but now this view seems doubtful to me because 
of the following contradiction. 

Indeed, in the referenced §9 of Begriffsschrift Frege had written, "On the other hand, it may be that the argument is 
determinate and the function indeterminate." This is idea can be explored in terms of predicates and sets, as Russell 
did. But it is also a fairly straightforward description of functional programming — the idea that functions should be 
treated as first-class mathematical objects (like numbers). This is implemented in several computer programming 
languages (most prominently Lisp |27|). Even beyond the conceptual appeal of such an object-oriented approach, 
it can be practical to solve systems of equations with functions as unknowns, when the type of function is suitably 
restricted (e.g. to polynomials over a given set of variables with coefficients in the finite Boolean field F2; Boole's 
function development method allows parametric descriptions of arbitrary logical functions of any desired arity). Some 
challenging problems in logic, such as SmuUyan's puzzles about liars and truth-tellers [251, amount to solving for 
unknown logical functions rather than unknown elementary truth values. The calculus of infinitesimals is another 
place in mathematics where functions are treated as first-class objects. What is integration except solving a system of 
constraints (on rates of change and perhaps boundary values) in which a function is the unknown variable? 
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C NOTES ON MODAL LOGIC: THINKING OUTSIDE THE BOX 



We can use dynamic polynomial logic to define modal operators in terms of solution-value sets. First the alethic 
modes corresponding to the members of 2^"- ' ^ (here using s to abbreviate the solution- value set (p) of the objective 

p{x\,X2,--- ,Xn) relative to the set A of axioms): 

o □ (or Di) necessarily true: 0{p) means i = {1} and D{—'p) means s = {0} 

o Dq necessarily false: □o(/') means s = {0} and □0(^/5) means s = {1} 

o txi ambiguous (mnemonic 00 for all truth values): cxi (p) or [xi (-ip) mean \s\ > 1 thus i = {0, 1} 

o unsatisfiable (mnemonic empty set 0): <Z>{p) or <Z>{-'p) mean i = {} 

Then some hybrids, each corresponding to several members of 2^*^'^}: 

o (or Oi) possibly true: (){p) means 1 G s thus s € {{1}, {0, 1}}; 
and <){-'p) means 0€s thus s € {{0}, {0, 1}} 

o Oo possibly false: ()q{p) means E s thus s S {{0}, {0, 1}}; 
and Ooi^p) means 1 e s thus s e {{1},{0, 1}} 

o □ definite (mnemonic 'get to the point'): □(p) or S{^p) mean \s\ = 1 thus s G {{0}, {1}} 

Perhaps we should use □(pjA) or 0(p\qi,q2, ■ ■ ■ ,q,n) etc. to make the axioms explicit. 

These operators and their negations cover 14 of the 16 members the power set 2^'" '' . Of the remaining 2 members 

{{}, {0}, { 1}, {0,1}} is tautological and {} is impossible (even the empty set is not a member of itself). Modal 
operations are recurrence relations; they must be evaluated as dynamical systems (though it is fine to report static 
interpretations). All these modal operators have the property that they can be evaluated during hypothetico-deductive 
analysis: given a computed solution-value set s it is possible to tell whether any operator is satisfied using standard 
set operations (equality, inequality, membership, complement, cardinality). For multiple variables (x\,X2, - ■ ■ ,Xn) we 
should probably use the solution set of all variables as the state of the dynamical system. 

For t/-valued logic we could introduce d distinct modal necessity operators Do, • • • , Dd-i with each meaning 
s = {k\. The remaining 2'' — t/ — 1 non-empty solution-value sets in the power set of the set of elementary values (each 
necessarily with cardinality > 1) would satisfy ixi (thus making it hybrid rather than singleton). We could likewise 
introduce d distinct modal possibility operators where each ()k means k€: s. Whatever special value k* is considered 
'true' would customarily be omitted from the operator sign: thus the unadorned □ and instead of n^* and Ojt* • 

Note that the position of negation matters. For example ^(D/?) means i 7^ {1} thus the solution-value set s could 
be {}, {0}, or {0, 1}; whereas \I\{-^p) means the solution-value set s is definitely {0}. In other words we distinguish 
between 'not necessarily true' and 'necessarily false'. Note also that 'possible' here does not follow classical modal 
logic: (){p) is not equivalent to ^□(^/?) because the latter includes the possibility s = {} that p is unsatisfiable whereas 
the former does not. That is, we distinguish between 'possible' and 'not necessarily false'. Also ixi {-'p) means the 
same as ix: [p): the negation of an ambiguous formula is itself ambiguous. 

In this framework □(/?) does not mean the same thing as p. In fact the axiom h □(/?) is quite a different assertion 
from the axiom h p; in isolation, the former specifies a dynamical system with a fixed point at the state that p is 
unsatisfiable. However the solution-value set ^(□(75)) evaluates to {1} exactly if the solution-value set also 
yields {1}; in this sense □(/?) is similar to p when these expressions are treated as objectives. 

Only unsatisfiability is 'explosive' : from a contradiction it follows that all formulas are unsatisfiable (including 
the constants and 1: e.g. ,i^{i=o}(l) := {1 : 1 = 0} ^ {}). Thus if (Z){p\A) is true for any formula p then it must 
be true for all formulas subject to the same axioms A. Also, the sense in which algebraic proof is monotonic is that 
adding constraints cannot add members to solution sets; it can leave them unchanged or remove members. Thus for 
any sets A and B of axioms and any formula p we have ,S^a\jb{p) ^ ^(/?) and .5^a\jb{p) ^ '^b{p)- It is not the case 
that □(p|A) guarantees D{p\A,B) since the additional axioms B could render the whole system unsatisfiable; thus we 
could have □(pjA) but (Z){p\A,B). Algebraic proof is global not local; all constraints must be taken into account, and 
the solution set y is what is left after all infeasible solution values have been ehminated. 
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D BY GEORGE! PROPOSED LOGIC QUERY LANGUAGE 



Sketch of commands for a proposed computer system that implements 'Logic Query Language' incorporating dynam- 
ical systems, Boolean translation from logic to algebra, and polynomial equations. 

/* george: computational algebraic logic (LQL) */ 
pareimeter c in {0,1,2}; 

b = 2; // simple macro assignment, resolved at compile time 

c := l$x|; // update rule: c gets size of solution-value set $x of x 

x"2 + b*x + c == 0; // equation template 

x"2 + b*x + |$x| == 0; // same constraint, implict parameter definition 
parameter t in FF(2); 

t := $y == {l}; // update rule: t gets true if y was proved a theorem 

t := ?y; // same update rule: t gets true if y was proved a theorem 

y == 1 - t; // template: y true if it was not proved a theorem 

|- y <-> !t; // same constraint in logical notation, |- for assert 

y == ($y != {1}); // same constraint with implicit parameter definition 

y == !?y; // same with modal operator ? (necessarily true) 

y = !?y; // interpret as static constraint or recurrence? 

// Fibonacci-like recurrence relation: update rules but no constraints 

parameter tl, t2 in NN; 

tl := t2; 

t2 := tl + t2; 

tl[0]=0; t2[0]=l; // initial conditions give tl[t] : (0,1,1,2,3,5,8,...) 

function F(2,2) indef inite<z> ; // arity 2, each argument in FF(2) 

// creates variables z[0,0], z[0,l], z[l,0], z[l,l] for coefficients 

// call as F(x,y), for which polynomial function of x, y, z[i,j] substituted 



paramieter d; 

function G(d,d,d) indef inite<w> ; // arity 3, each argument in FF(d) 
// creates variables w[0] . . w[d"3] for coefficients 



fimction G(x,y) : x -> y | | y -> x; //a definite function 
G(x,y) : X -> y I I y -> x; //a definite function 

(:A x,y,z: F(x,y) -> z) // imiversal queuitifier 
( :E X : F(x,x) -> 1 ) // existential quantifier 

parameter 'in in FF(2); // keyword as identifier 



7o solve x; // runtime query solution- value set for x 

% $x; // same: returns recurrence F(x), domain U={0,1,2}, objective x 

% ?x; // ask if dynamic theorem: exactly one fixed pt x == F(x) with x == 1? 

7. $ tl (0,1); // to get sequence (0,1,1,2,3,5,8,...) 

% $ tl; // to get collection { ■C0}@(0,0) } 

7. $ X @@ // collection of sequences $x@0, $x@l, $x@2 
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E WHAT THE TORTOISE SAID TO ACHILLES ABOUT WHAT THE TORTOISE SAID TO ACHILLES 



Achilles had once again overtaken the tortoise, and had seated himself comfortably on its back. 

"So you've understood the Liar?" said the Tortoise, "even though he told you he was lying? I thought some 
mathematician or other proved that the thing couldn't be done?" 

"It can be done," said Achilles, "It has been done. Quod erat faciendum. You see, it was a simple matter of solving 
a system of simultaneous equations; and so — " 

"But what if the equations did not happen all at the same time?" the Tortoise interrupted, "What then?" 

"Then I shouldn't have solved them," Achilles modestly replied; "and you would have got several times round the 
world, by this time!" 

"You impress me — compress, 1 mean," said the Tortoise; "for you are truly dense, and no mistake! Well now, 
would you like to hear of some logical propositions that are supposed to be true or false, while they are really neither 
of the two?" 

"Very much indeed!" said the Grecian warrior, whereupon he produced a notebook, a fine quill pen, and a small 
pot of ink from his shield-pocket protector. (He wondered to himself, "Am / not the protector of this pocket?"). 
"Now listen carefully, for I shall reveal to you the Principle of the Excluded Muddle." 

Achilles rose to his feet and balanced atop the Tortoise, still writing furiously. "So my answer depends upon itself, 
rather Uke a turtle standing on its own shell?" 

"But just what do you suppose that answer depends upon, my gallant Greek? Surely you can tell me something 
more particular about it." 

"Very clever, young Tortoise, very clever" said Achilles, "but it's turtles all the way down! I dare say any particular 
answer would be nonsense." 

"Nonsense indeed," replied the Tortoise, "But which kind of nonsense? There are two different methods to this 
madness, the unsteady kind and the contingent kind. You can tell them apart if you start from the beginning, which is 
after all why we call it so." 

"Let us start with the very first turtle — or Mrtle if you like, since it has no beginning." 

Achilles sighed, "With such a long and wandering explanation, which will surely vex some Logicians of the 
Twenty-First Century — would you mind renaming yourself Tor-tu-ousT' 

"As you please!" replied the rational reptile, admiring the warrior's tidy notes, "Provided that you, for your 
excellent permianship, will call yourself A Quill Easel" 
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